


f 


ns Cain 






Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1986-09 


Texas Instruments 4100 GPS positioning software 


Brown, Kevin T. 


Monterey, California. Naval Postgraduate School 


http://ndl.handle.net/10945/22146 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
D U DLEY research materials and institutional publications created by the NPS community. 
get Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS'‘s first 
KNOX appointed — and published — scholarly author. 


LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


Pe ts 


Lamsew 0 vane oy 


Boe ds bat tae 
eA al ana , 
De tude? san: 49,07, 
Sao G2» wee Fils 
SR del eng gas ees 


rr ee ee oar) ay 
rey , 


a ha Fae er 
marta, 


Pieet eas 


to Wag y trike 
Ttmirs ge 


Ne hae es Ng 

a Re. 2%. Vaart ad’, 
“ae EF oeagis! SOD BF 6 a2 
6 stage oon, 

ae fo e%e 


La 

2) 84 bth 
Lard So MEY 
al font th 


Lhe] 
OF ar are # Arm tears 
Nw tet Chae a” 
“hoy ty 
nate 


Viera eadeneg, 
Reatageg Sint 
Mera 
we ate 


EERE easy 
em stote* an sad ye fates 
AR? 6M to 


Fin! i a= rf fale Pad 
OW he 

Cie or beaserg tle 
Ata. ec@ ASR 


154 Ae, 


Thy CT 


a ee 


‘= Mi Acrat | ue Ae tian 
iS th te re, 34 detcashiags r te 

fk, ahaha: ca tthe chsaenhegn” 

JP bintsetrdcady 
As Foe g¢ 
45 15 603 vedas 
dashes Vaasa 
<M” bias =A Sees, 
inwea Phe 

* 0° Rihcoe 


e 
2 aes 


Le i ‘a. 
£3 Pe ee Fy 8 hentia he wa, cal t26'O-® 
SiR tha ae sas ges 
Wee be 
Rd gar ob 


ro 
tos his a vat 


re 
TS Ole gs 
fh Fa SASGAZs ot 
Treats dae e: 4 
«4%. Vita trang 6b! te 
Lea ee ae Pe 


UGA ey 
wwiye? 


BF Fat Oe og 
RA tel tot tig da 28% ang 
SUR sy fag MEI NALA ges, 

2 obi” Omaha 5 Ow At hal, 

45 ER Casha Ory g? 


2" ita VAes 
WTR sds ac bate 
rf & Orn! o- deep 
tai tens 
mT et ot te 


Yee ase ote 

*e Gomi ed gat 

OFLA INA none 
Viet te aderaca 


shaban 
bteta,o* 


Nitetng yg 
alenns a 


Weta ag, 
th args 
: 


Leg, 
Pan 
eh 


ean eek 
| MlEGs 


pleins 8S ace 
*EA LASTER hog 
han 


te aiesd ES 
Shera MS orw b 


Rig toa 2. eS.) i 
Ms ale. oar) rie, 
4 Phage 

©4460. a4: ,, 


20 5-4 alas Mzha 
ieee; 


ete oo) x 
tg nach § MoM f 
ea z % € 4 
07h 6 


Atenas o¢ ' A 


S 
*, 0 alece: 
Sh Arete MG a at, 3 
Whee iy 218A on, 
‘maa 8 ate Om 
Are a2 Ssake—acbran han tol 
Me BS Shee *s ie case 
WiAl 4s 82 Luan ca cdeene oe tase 


a og 
A thetr ican 


a; dea. 
© ona 
*asere- ar 
te: he Ris deg ie 
edad 


S En am z 
Mase 
ee 
Singh ¢ 
fasice. 
eles d Rody 89 ote. by: 


ing 


% 

- WR itarern! 

Cand Sa WEA tog 
Syedet he tars: 

beet 


Wey WES gt SsBae® tee 
es "ens 


hath ts ed 
easels . 


eee 
oera np 


Wie tet ne 
Oem OR ian eden hy Aste as 

. APS cated, a2 

ww rr * j S 

? sag 

Te he 


Wes hekng 
oe ge Si WEw ly 
ork © ote 


La 

Rwy ad aw 

1 ‘ i? s ™ 

f res’ fd pee 

aara, ? 

5 Soy eee Vat aA s ; 

* . * ise as x 

Niel 5 aD | } 

: 7s i } 

Hote, Pah weurttirs 

$+ APN Paes wren U4 Z 
; 


‘ Paes) ert | 
rer: bg 
a ery t ded Uni Mews sa ceanareege’s, A 
wegt, Lay cits ‘ 1 
ahd wtih 5 te A oeer 


eye 


re ary Yavatuyn Mie, wy 0 
wForytie. Potter arn ey % 
2 Woes 
ee ee es Mare Yi 
Ftd et yt yeg er” =! ve 
BOUT Sere “ATW Oe ance 
thease fy 


BENT as Be eeitat 
atte dc : " 
elaine 

cela yy! 
BOL vag 4 bd 


ch an Fap (oTin 
reek tae ih -c¥la « 
ae aa lry e 
ets tow on ad. 

Ld Vews timing: 
ETL Le 2 
td RSs, eg, whe 
Ve atee ce eid esag 
. 
oe *y indi 


aragent 
WN seen! 


SAF oe, 
‘ fs tenes 
era? Fe * 
“data og ; 
2 S39 “tht ate ye, 
Faye AF eta as 
hay releoyt 
Sa atar Fo due’, 
ycSetee 
i} 


ba 
She ca loe. Mate 
A fade raet, 


wth aatagl 
Peet ae ee ea Y 
te lt sme 
Pier e berg te! 


bo eed 
Me i 

- 4 ead vend, 

LT er Tey 24 Rete ata tes clara. 

Te ate = Hn gh Tees wartisge 

Gr ae Notas sett ae 

ose 1s None oy 


ke ee 
Shin dh Cedar ing 
=f cao 


ear tie che Vets ater tay 
We Ver ae) Patt Caraws *% 
grinds yea, 4 ston tmaeic $< an wg 
k Rave d nists 
ar 
THAMER eye ira . 
Ca dee Se - i & me) 
te reaedeaas fy adarle cy? ue Nears at 3 7 oe Pe ile Le 
dat heiecninas rH ‘ . .: Atay FEE i ee 3 99 Le Chath hh Seal 
‘ EY a s4y Fein tat ' = > Se td tet y is at 7 ¢ Pee he ueks ace ae : SS ets 
PuReeLL Se Tero ga ae ; = v8 asada ie Se PD 2%? we nd as g FS oy as itig 4 5 r 
ERG Oe ope 5 é We Tbe re eary oF ‘ ELA An ag 4 f i rigs : 7 ‘ 2 . 
state BPS entnsnt ls Bay fate’ dad 2 tyne j 
fsa “ : eh ie 
* * * 
= a dta ; 
i $9 kz HG ay* 
Ps eetyn Cg Ie 
Lay & fetuea 


. 
3 Lhe 
onto ob Bex Ht 
eM ecw ala Seasotag 
windadary eaehe 
Frat tianw Papa 
Stes TL bg 


Tatas 
we dees Fics 
em Tr wheats 
whey 


afte ae 


as 


ei tn 
: Pi Beta Care tr 
i 7 4 
ae hey ae 
seigte Secst ‘ 
+2 fae Meigs ay 
Le eae b Sai 
2°? eieey me ek 
Atphene t sigierstes 
Sot hte tg ee 
te Td hetes 
ge2a 


tueie 
otts}, 


Pega e 
£80 %s gp, Mgxe 
MALICE tee F 
irs TRS capes 
ede Sadbee ae de et wz 
TM Rewtee eg, 
ae segs sai). Ra. ite 
Shek sts ate GN So 
tetserte Ht eos 2 Rgt et ges 
t C=. a8 
thd Ged” Hal ae eet Rs ate 
wr tt sin 


=e ies 
tae coon ce% 
Pee 


sy 


a 
nt aS oe id ‘| SS 
SHEEP te, oe 

Hety! Be nd 

ptt dtae pies bat 

hl ete Upalt 


£* 2° 


7h. ee 


pi tees seas gy’ 
* tf R i 


NAS MRSS Pe 
SHS Bet ate ces 
af. 
‘ ri : 
Per, f 
; ‘ , 4 Y, at eC-4T pe: 
“ c i . 
Kar a : be A ‘ 
. hoy ‘ e . teu 
Pe ae oo) 
: 


# 


=e 
ty ens 


am 
3h 


Mtn ut, + 4 . . qt 
Me ete Gt . Chiesa ed, 3 at 
ay a" é é a Zi - < 


‘S” SS aK es ghee, 
fe re Ft 


vias 
i Oh 
J 0s tas " 


2 aide 
Ph 


i 
ssa 
« 


ro 
ay rs 


¢ 
3" 2"e"iges 
ig PMS fpr gg I 
° saree afekeng 
: what Ey 
; Vere ely ; 
Fe ate . agen et"p> 5&4 A * my 
re aE a ae yal tt 
SEU bs Ie ad a ates feeke’y He" Cin nrg racy 
MEPS adh eth Rea if ve * Ag" varsat 
Paty anne gta BY gt 5” .2° t}¢< Gj Pra se ae 
EO Ar ae 2 et ee 


Pel pure acy 
ARE alate 
an oe 
"Pure 
Efe ay 


797 49 5% 
’ ars 


we 
: tae PEP ae 
LEP gre la: bthest Beret ho] 
la a 
ae ag 


~ 4: 142° drags, 


ry So eee 
are l erg a se 
Se Egg: 
ser 


WV tee 


“a ghasy 

ee feiss a3 ee Lae ey * 
‘Fe SBC Noe ee \ 

7 Pons abegt eeege ge O73 elwtat “¥r3 a) Pith 

> x . # te pyar g ' Senate Weer ey if 
the py “+ Dye wey ° iF Wa sha, 

AR CL a ta K Eee 

es J big he 


°F MN agy * 


a 


F atirepanta? 
ee Pecans faty sty 

Dy a A rs ee = gh gees 
“ia soy ede iy ba bo netstat 
ye ye 


f etareg? 
oil eey 

“st , 
oct 608 a 


. a 


: 3 see 
ye. “" < SOE oD rep are Ly we 
He Learn , * prea? y f 9? 3 "te ited VESTA o-0i areeprere gs whales 
Pad Mak rie tu: i pis ‘ A SF < Nays Pwermetare rte’ Tere 
fe i) ee dad Wipe, : 4°¢ © Samy ~yt ores os 
4 é Dil dei 7 WES ss CS 
ms f dere ef ? 2 =r 
Fs try wtery eet ycy! Setar ie rer eres " 
STASI Han eRe or att WSs saaielat F,Weer9 aa geet tohy 
P17 sprays : xs Ms fy: es z hee 
new ry 


sy 
ert ie 


7 & i 
US pes te, 
OR SN | 
ra ie 4 ‘ie 
"fa pn'g + 
asya'y 


ety 
Loar Mar SPH 


es 
Pe Yk 


re 
i] 


er 
4 Eat 


} : a AA ones Ty 
iP "¢ pyy oe pea ‘ “ae f : A tel ue x va 
sxye hide Sey Ee ; oan utes wet ratio tete 
Wh die de A ot rey el Ap “9 
Raa 4° Ad Gas rae oo Re Las 

ee WN ten, ¥ n, 2>3"o 





DUDLEY st02 LiRRARY 
WAVAL POS" 7t4 20 sei SCHOOL 
MONTERSY Grou mene 95945-5008 








NAVAL POSTGRADUATE StHOOL 


Monterey, Galifornia 





THESIS 


TEXAS INSTRUMENTS 4100 GPS 
POSITIONING SOFTWARE 





by 


Kevin T. Brown 


September 1986 


Thesis Advisor: Muneendra Kumar 


Peeme ved Tor public releéase- distribution is unlimited. 


1230148 





SECURITY CLASSIFICATION OF TH.5 PAGE 
REPORT 


Ya REPORT SECURITY CLASSIFICATION 


unclassified 
2a SECURITY CLASSIFICATION AUTHORITY 


2b DECLASSIFICATION / DOWNGRADING SCHEDULE 


& PERFORMING ORGANIZATION REPORT NUMBER(S) 


6a. NAME OF PERFORMING ORGANIZATION 





Naval Postgraduate School code 





6c ADORESS (City, State, and ZIP Code) 


Monterey, California 93943-5000 





NAME OF FUNDING/ SPONSORING 
ORGANIZATION 


8a 


8c. ADDRESS (City, State, ard ZIP Code) 


11 TITLE (include Security Classification) 


TEXAS INSTRUMENTS 


6b OFFICE SYMBOL 
(if applicable) 





8b. OFFICE SYMBOL 
(lf applicable) 


DOCUMENTATION PAGE 


to. RESTRICTIVE MARKINGS 


3 DISTRIBUTION / AVAILABILITY OF REPORT 
Approved for public release; 
unlimited 


distribution 


5 MONITORING ORGANIZATION REPORT NUMBER(S) 


7a NAME OF MONITORING ORGANIZATION 






Naval Postgraduate School 
7b. ADDRESS (City, State, and ZIP Code) 
Monterey, California 93943-5000 


9 PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 







10 SOURCE OF FUNDING NUMBERS 


PROGRAM 
ELEMENT NO 


TASK 
NO 


WORK JNI™ 
ACCESSION NC 


PROJECT 
NO 





4100 GPS POSITIONING SOFTWARE (UNCLASSIFIED) 





*2 PERSONAL AUTHOR(S) 


fli 


Brown, Kevin 











135 TIME COVERED 
FROM 


3a TYFE OF REPORT 
Master's Thesis 





TO 








"6 SLFPLEMENTARY NOTATION 


COSATI CODES 


: A 






1§ SAGE COUNT 
8 0 






14 DATE OF REPORT (Year, Month, Day) 
1986 September 


18 SUBJECT TERMS (Continue on reverse if necessary and identify by block number) 


Global Positioning System, Texas Instruments 4100 GEOST 
a Receiver, CDC Cyber 865 Computer, IBM 3033 Computer, 
foe ee aco catic ~~ Dynamic Positioning 


"9 ABSTRAC™ (Continue on reverse if necessary and identify by block number) 


The Hydrographic Sciences Group of the Naval Postgraduate School's Department of 

Oceanography has a goal to establish and continue a research base related to the 

Global Positioning System (GPS). To meet this goal a library of programs capable 
of reducing GPS satellite data to position and time solutions, independent of the 
receiver's navigation solution, is being acquired. 

This thesis describes the modification, documentation, and establishment of Texas 
Instruments 4100 GEOSTAR GPS software written by the Naval Surface Weapons Center 
mer their CD@ Cyber 665 computer,~so that it will run on the IBM 3033 computer at 
NPS. Position and time solutions for either a static or dynamic receiver can be 

achieved using this software. 


21 ABSTRACT SECURITY CLASSIFICATION 
unclassified 

22D TELEPHONE (Include Area Code) 

(301) 227-2211 

83 APR edition may be used until exhausted 
Aflotner editions are oosolete 


i 


7S @3UTION/AVAILABILITY OF ABSTRACT 
RZ) ONCLASSIFIEDCUUNUIMITEO (1 SAME aS 2PT 
C24 SAE OF RESPONSIBLE |NOIVIOUAL 
Dr. Muneendra KAumar 


OD FORM 1473, a4 Maz 


4 
pa Ne 


Moric users | 











22¢ OFFICE 


GS 


SYM3O. 


SECURITY CLASSIFICATION OF "41S PAGE _ 


Approved for public release; distribution 1s unlimited. 


Texas Instruments 4100 GPS 
Positioning Software 


by 


Kevin T. Brown 
Cartographer, Defensé Mapping Agency 
B.A., University of North Carolina, 1982 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN HYDROGRAPHIC SCIENCES 


from the 


NAVAL POSTGRADUATE SCHOOL 
September 1986 


ABSTRACT 


The Hydrographic Sciences Group of the Naval Postgraduate School's 
Department of Oceanography has a goal to establish and continue a research base 
related to the Global Positioning System (GPS). To meet this goal a library of 
programs capable of reducing GPS satellite data to position and time solutions, 
independent of the receiver’s navigation solution, is being acquired. 

This thesis describes the modification, documentation, and establishment of 
Texas Instruments 4100 GEOSTAR GPS software written by the Naval Surface 
Weapons Center for their CDC Cyber 865 computer, so that it will run on the IBM 
3033 computer at NPS. Position and time solutions for either a static or dynamic 


receiver can be achieved using this software. 


Lite 


TABLE OF *GORWE NTS 


INTRODUCTION «o.oo sn copes sae 9 
BACKGROUND 6 siseosw 5 Ree eo 0 reese nese 10 
A. GLOBAL POSITIONING SYS TEM eee 10 
I, Space Segment .<.2 . 2. :c0ge eee 0 occa ee tack 10 
2. Control Segnvenmt [iee. =. eres a ee nr Il 
3. User Segment <2 0.05.% .. «0. aeepeeenee te cee emer cr I] 
B. TEXAS INSTRUMENTS 4100 NAVSTAR RECEIVER 
(TT-41OO)> 5 is ves ce ce Bae eo, oe eee ee igi 
I. Characteristics and Descriyptroi a ieeeee tee eee nner reer I] 
2. Receiver Measurements and Postion Calculations... 12 
3. Position Accuracy «2.2 2 fevers ycecec ls cee tee eee ie 
C PSEUDORANGE DATA AND THEIR 
CHARACTERIS@1G@siae ees ee 12 
D. STANDARDIZED EXCHANGE FORMAT FOR 
NAVSTAR GPS GEODETIC OA 3. ee eee 16 
I. The Basic Pormat ........5.4. 5. ele eee eee 16 
2. Coordinate SyStenny 5c 5c spe eesnern ere reece 17 
DATA PROCESSING aac. « «20s. 5 < See [8 
A. PRE-PROCESSING }.......2..:. seers re [8 
I. Manipulation of Cassette Recorded ara... eee ee I8 
2. Alignment of Common Time-Vaecceds eee 18 
B. DETERMINATION OF PSEUDORANGE POSITION 
10) FO ip (O'R 6 sO Comb eo boo boom ous 22 
I. Static Positioniige occ... ia eat eter pee teks ae we 
2. Dynamic Positrons 50.2 ores: eeee meee ere eee tae eer naa 34 
TEST RUN RESULTS AND VAI Adler ee 43 
A. STATIC POSITION SOLU WO Nese tele eee eee ee 43 
B. DYNAMIC POSITION SOLUTION Ribs Gitte rere 43 


Ne Sh CERNE oc! a6 nS nn 46 
VI. Eee OVI IZ sO RE. sk a eee ee ee eee 47 


APPENDIX A: CONSTR SUBROUTINE DESCRIPTIONS AND 
SES) UN CS SR a coal) Coe ee 48 


ae Neel nee lyre rOk PROGRAM CONTR 2... ees 50 


per eC: TPE XAS INSTRUMENTS GPS SATELLITE 
Rese) We TPE) 2 JE OOROAY 8 00 Dae a rn = 
APPENDIX D: enn SUBROUTINE DESCRIPTIONS AND 


we OoUN ED UNPUT ee HN SS TEES POR 


APPENDIX F: Se ee SUBROUTINE DESCRIPTIONS AND 


INSIEISS LUNG? 26 Son or geen tore teers a le er 66 
ari G: TNPUT GUIDE FOR PROGRAM KALMNZ2 ...........6... 69 
Proreerweoi x Ti ABBREVIATIONS AND ACRONYMS ........6..0c000eee- us 
Lowe | COLE [EJ STEIRIE INC |p) SU St har en heer arr a re 76 
Pee ee PAA LOIN EIS eee ee ee eee. es st 2 ii 


ITT. 


LISIWOR Abies 


Tl-4100 POSITION ACCC NG a ee 


CONITR AND SOLTOM SOLUTION RESULTS 
CONSTR AND KALMN?2 SOLUTION RESULTS 


* «© «© @# « ¢ @ @ @ © @ @ © @ © &© 8 © 


oe ese 868 8 @ 8 © @ @ © @© © © # @ 8 8 


e*eeeeee#eeeeee#ee*e# ¢e¢s8 e# 


21 


Sell 


eo OPP iGURES 


TI-4100 NAVSTAR System (From Texas Instruments, TI-4100 
CONE SAMMI ce te tc w esse vee eeaee 


blige tdntesor Processing 11-4100 GEOSTAR GPS data ............26. 


ACKNOWLEDGEMENTS 


I would like to express my gratitude and thanks to Stan Meyerhoff of the Naval 
Surface Weapons Center, whose help, by telephone and in person gave me the tools to 
make this thesis possible. I would also like to thank my wife, Beth and daughter, Anna 
whose love and support encouraged me to complete this work. Most of all I 


acknowledge my LORD, who gives all knowledge and wisdom. 


1. INTRODUCTION 


The Hydrographic Sciences Group at the Naval Postgraduate School (NPS), 
conducted its first Sea Floor Benchmark Experiment in May 1985. The experiment 
basically was to deploy a benchmark on the ocean bottom and to establish its position 
with geodetic accuracy. This was to be accomplished using GPS and acoustic range 
data (Kumar, et al., 1984). 

Bevond the goal of the Sea Floor Benchmark Experiment in May 1985 was a 
higher goal of establishing and continuing a GPS research base at NPS. In order to 
reach this goal it has become evident that a hbrary of GPS prograins should be 
established at NPS capable of reducing GPS data, independent of outside sources, into 
a format compatible for analvsis. 

The Texas Instruments receiver (TI-4100) was selected as one of the GPS 
positioning receivers to acquire data for the experiment. One version of the T1-4100 1s 
for commercial users and the other (TI-4100 GEOSTAR) 1s a tri-agency version for use 
by the Defense Mapping Agency (DMA), National Oceanic and Atmospheric 
Administration (NOAA), and the United States Geological Survey (USGS). TI-4100 
GEOSTAR software has been obtained for use at NPS from NOAA and the Naval 
Surface Weapons Center (NSWC) in Dahlgren, Virginia, for compatbility 
modifications with NPS’s IBM 3033 mainframe computer. Of this software, three 
programs from NSWC, written chiefly by S$. L. Mcyerhoff for their CDC Cvyber-865 
computer, reformat GPS data and then compute receiver positions and time bv either 
an iterative least squares approach or a Kalman filter. 

The main objectives of this thesis are to describe the effort in modifving NSWC’s 
software codes and algorithms for compatibility, and establishing them on the 
IBM-3033 mainframe computer at NPS. In addition, discussions of the mathematical 
models, results of test runs and validation of the modified software have been included. 
This software was acquired from NSWC for exclusive use at NPS. Future users at 


NPS, are encouraged to test and verify their solutions as correct before accepting them. 


Il. BACKGROUND 


A. GLOBAL POSITIONING SYSTEM 

GPS is a universal satellite positioning system that is independent of weather, 
time, or geographical position. It provides both position and time in an earth-fixed 
geocentric cartesian coordinate system as well as velocity data. Since it is a space 
based radio navigation system, it can provide accurate information to an unlimited 
number of users globally on or near the earth’s surface. Three-dimensional position 
data to within a 16-m spherical error probability (SEP), velocity to within 0.06-0.15 m/s 
and GPS system time to within 25 ns are possible instantaneously to properly equipped 
users anywhere within 500 nmi of the earth (Milliken and Zoller, 1980). 

The GPS is comprised of three major segments, the Space, the Control, and the 
User 

1. Space Segment 

The GPS when fully operational will consist of 18 active satellites (Space 
Vehicles or SVs) in six planes with near circular orbits at altitudes of 10,900 nmi and 
periods of 12 hours. Each of the six orbital planes is inclined 55 degrees to the equator 
and is to contain three SVs equally spaced. The configuration of the SVs in the orbital 
planes is such that a minimum of four SVs will be available to the user at all times, 
anywhere on or near the earth’s surface. Each individual SV transmits its broadcast 
ephemeris and clock correction, plus an almanac which contains orbital parameters and 
clock correction estimates for all other SVs in the system. 

Each SV transmits two radio frequencies simultaneously, L1 and L2 
respectively, containing navigation information in the navigation data message. The 
Ll frequency is centered on 1575.42 MHz and the L2 frequency is centered on 1227.60 
MHz. These signals are in a frequency range that provides good all weather operation 
and allows the user to determine ionospheric propagation delay. The L1 signal is 
modulated by a 10.23 MIIz clock rate Precision (P) code and a 1.023 MIIz 
Course/Acquisition (C/A) code (Spilker, 1980). The P-code is a random sequence 
binary code that provides the user the capability for high-precision positioning and is 
highly resistant to electronic and multipath interference. The C/A-code gives all users 


rapid acquisition capability and acts as an aid to gain access to the P-code. The 


10 


navigation data message allows the user’s receiver to calculate position, velocity and 
GPS time by providing broadcast ephemeris, clock corrections and almanacs for all 
SVs. The L2 signal is modulated by the P-code only. 
2. Control Segment 

The control segment consists of a master control station, a number of monitor 
stations and ground antennae located around the world. The USAF master control 
station is at Colorado Springs, Co., and monitor stations are located on Kwajalein 
Atoll, Ascension Island, Hawai and Diego Garcia. DMA also maintains three stations 
in the United Kingdom, Australia and Argentina, and stations are proposed for 
Ecuador and Bahrain. The imonitor stations acquire ranging data for each SV in a 
passive mode. These data are then relayed to the master control station where they 
are processed to determine the SV’s position and signal data accuracy (Wooden, 1985). 
The master control station takes the relayed data, computes errors, generates a new 
navigation data message and then uploads this message to each SV via the ground 
antennae. The ground antennae transmit as well as receive satellite control 
information. 

3. User Segment 

This segment consists of stationary, low, medium and high dynamic receivers 
designed with specific requirements for individual users. The user receivers are 
designed to receive and process SV data from four SVs either simultaneously or 
sequentially. The TI-4100 receiver determines the range from the SV by measuring the 
time delay of the SV’s specific epoch and the user generated codes. Then making 
separate ranging and carrier frequency Doppler shift measurements to four SVs, the 
receiver calculates the user’s position, clock bias error, clock drift and velocity. The 
position coordinates are in the World Geodetic System (WGS) 1972, an earth-fixed 


geocentric cartesian coordinate systein. 


B. TEXAS INSTRUMENTS 4100 NAVSTAR RECEIVER (TI-4100) 
1. Characteristics and Description 
The TI-4100 receiver is portable for field operations (it is about the size of a 
small suitcase) but can be mounted in a standard 19-inch rack if desired. A small 
lightweight hand-held Control Display Unit, contains a keyboard and a display 
window. The receiver is modular in design and allows for easy removal and 


replacement of circuit boards. 


1] 


The T1-4100 is a single channel dual-frequency, multiplexing NAVSTAR 
receiver and navigation processor that tracks up to four SVs (Texas Instruments, 
1983). The dual frequency receiver allows for tracking LI and L2 frequencies for cach 
of the four SVs. It also has an optional dual cassette recorder for recording GPS/SV 
data. If GPS/SV data is to be acquired in a degraded mode (i.e. less than four SVs) 
there is the capability of adding an external atomic clock time standard input to the 
T1-4100 receiver. The T1-4100 receiver, Control Display Unit, antenna and cassette 
recorder are shown in Figure 2.1. 

2. Receiver Measurements and Postion Calculations 

The basic measurements made by the TI-4100 receiver are the arrival times of 
the SV-generated signals compared to the receiver clock, which contains a stable 
oscillator and a counter (Texas Instruments, 1983). Four pseudorange measurements 
plus the navigation data message provide enough information to calculate a 
three-dimensiona! user position solution plus the user clock correction to true GPS 
time. The solutions are obtained using four observation equations and four unknowns. 
Doppler shifts due to the relative velocity between the SV and the user are measured. 
Measurements are also made of the phase differences for the LI and L2 signals 
between the SV and the user. | 

3. Position Accuracy 

Error statistics for the TI-4100 receiver are given in Table I. It should be 
noted for this table that Operating Modes using time aiding, have a positive growth 
rate error of 2.0 to 10.0 m/h (l-sigma) for frequency standards that allow navigation 
calculations to estimate precise frequency (Texas Instruments, 1983). The positioning 
accuracies in Table I were based on the TI-4100 receiver observing pscudoranges, 
Doppler shift data, phase differences for LI and L2, then computing position, time and 
velocity solutions. The pseudoranges from the Tl-4100 receivers are the major 


observables considered in this thesis. 


C. PSEUDORANGE DATA AND TITEIR Chek ACTE apices 

The basic operation of the GPS depends on the user determining pseudorange 
and range rate to a number of GPS SVs which have precisely known ephemerides 
(Texas Instruments, 1983). Pscudorange is the apparent range, which includes the 
clock errors between the SV and the user’s receiver. The transit time of the navigation 


signal between each SV and the user is measured by the receiver and then scaled by the 
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_ Figure 2.1 T1-4100 NAVSTAR System 
(From Texas Instruments, T1-4100 Owner’s Manual, 1983). 
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TABLE | 
T1-4100 POSITION ACCURACY 


(From Texas Instruments, TI-4100 Owner’s Manual, 1983) 


1-Sigma 1-Sigma 
Position Error Velocity 
Emnot. (m) Ero 
Operating Mode Stawist les p C/A (m/s 


Four satellites 


Stationary 

0.5 g acceleration 

2 g acceleration 

4 q acceleration 
Ditferential 
Stationary differential 


WWWWWL) 
J 


1 
OVOVOODOO 


Three satellites 
Degraded (altitude aiding) Horizontal 


Differential degraded Horizontal 
altitude aiding) 


Degraded (external cesium 350 
frequency standard time 
aiding 


Differential Degraded | 221) 
external cesium frequency 
Standard time aiding 


Two Satellites 


Degraded (altitude and Horizontal 
time aiding) 


Differential degraded Horizontal 
altitude and time aiding) 





speed of light to compute a raw pseudorange. Each raw pseudorange must be 
corrected for tropospheric and ionospheric propagation delays, SV clock offset, user 
clock bias, relativity corrections, and earth rotation corrections. The pseudorange 
calculation 1s given by Equation (2.1). 

ne = KS =: cAta; ae c(At,, - Atc:) (2.1) 


where i = 1 to 4 (pertaining to four observed SV’s) 
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R;* = pseudorange to the SV 

R; = true range 

c = the speed of light 

Atc; = SV; clock offset from GPS system time 
At,, = user clock offset from GPS system time 
Atay; = propagation delays and other errors 


The tropospheric delay is estimated using simple refraction models with standard 
atmospheric values. Ionospheric delay is calculated by taking advantage of the dual 
frequency measurements of the LI and L2 signals. The SV clock error (difference 
between SV clock and GPS system time) ts modeled by the Control Segment and these 
correction coefficients are transmitted in the navigation data message from the SV to 
the user. User clock bias is due to the fact that the receiver clock is not perfectly 
synchronized to the SV clock. Relativity corrections are due to the eccentricity of the 
SV’s orbit. Earth rotation corrections are due to the motion of the earth during the 
signal’s propagation from SV to the user. Once four raw pseudoranges have been 
Comected “or esystematremm@errors, etheyeware *usedeetorcaleulate «the» receiver's 
three-dimensional position and time. 

The precision of position calculated by pseudoranges is dependent on the 
geometry of the SV’s with respect to the user. This geometry is continually changing, 
even if the user’s receiver is not in motion relative to the earth. If the intersection of 
the pseudoranges 1s orthogonal, errors in the pscudorange measurement are minimized 
whereas if the angle of intersection 1s small, errors are large because of the geometric 
relationship. 

Geometric Dilution of Precision (GDOP) is a composite “indicator” expressing 
the influence of the SV’s geometry on the accuracy of the user’s clock offset and 
position determinations (Milliken and Zoller, 1980). Other terms to describe dilution 
of precision in position and time are: Position Dilution of Precision (PDOP) which 
yields dilution of precision in three-dimensions; Horizontal Dilution of Precision 
(H}DOP) which gives dilution of precision in the two horizontal dimensions; Vertical 
Dilution of Precision (VDOP) which yields the dilution of precision in the vertical 
dimension; and finally Time Dilution of Precision (TDOP) is the dilution of precision in 
the user’s tine bias (Milliken and Zoller, 1980). GDOP, which is dimensionless, 1s 
calculated by I:quation (2.2), (Milliken and Zoller, 1980), 


GDOP = ((PDOP)* + (TDOP)?|/4 (2.2) 


The best geometry for a four SV configuration is the combination which results in the 
lowest GDOP. 

Once pseudorange, Doppler shift and phase data have been acquired, and the 
receiver's field solutions calculated, post-processing of these data may yield improved 
position and time solutions. To make TI-4100 data accessible to other users with 
different computer processing systems, it has become necessary to format SV data in a 
manner that can be exchanged among independent users. 

D. STANDARDIZED EXCHANGE FORMAT FOR NAVSTAR- GPS 

GEODETIC DATA 

This Standardized Exchange Format (SEF) was proposed by V. Dan Scott and J. 
Gary Peters (1983), and was adopted to make possible the universal exchange of 
NAVSTAR GPS data. To do this, serious compromises were made in regard to 
optimum magnetic tape utilization, computer execution efficiency, and computer 
program size (Scott and Peters, 1983). These compromises were required to allow for 
maximum tape compatibilty between different computer systems. 

1. The Basic Format 

The building block for each standardized data exchange tape ts the record, 
which 1s 80 characters. Under this format, an unlabeled 9-track, ASCII-coded tape 
with an 80-character blocking factor 1s used. The header (or the first) file on a tape 
consists of comment records and is not considered part of the data sect but 1s used to 
help the user to identify the contents of the tape. When the last data file is written on 
the tape a special end of information file is written (Scott and Peters, 1983). 

On the tape, records fall into two categories: Control Information or Data 
Information records. Control Information records are used to identify which 
FORTRAN Format statement to use when reading data records. Data Information 
records consist of one of eight data items that are uniquely identified by a four digit 


data item number. The eight data items under the SEF are listed below: 


(1) Constant equipment and software specific data 
(2) Campaign and station specific data 

(3) Calibration data 

(4) Time tag update data 

(5) Space Vehicle specific data 
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(6) Weather data 
(7) Real time solution data 
(8) Measurement data 
The goal of this format is to be independent of hardware and software design. In 
addition to data formats, the reference Coordinate System for the data must be defined. 
2. Coordinate System 
The WGS 72 an earth-fixed geocentric cartesian coordinate system is defined 


as follows (Seppelin, 1974): 
X-Axis = Intersection of the WGS 72 reference meridian plane 


and the plane of the Equator. 


Y-Axis = Measured in the plane of the Equator 90 degrees 


east of the X-AxiS. 


Z-Axis = Parallel to the direction of the conventional 


international origin (CIO) for polar motion. 


The reference ellipsoid for WGS 72 is defined by the parameters semi-major axis, a = 


6,378,155 m, and flattening, f = 1/298.26 (Seppelin, 1974). 


HW. DATA PROCESSING 


The pre-processing of GPS data is carried out to obtain a format suitable for 
further manipulation. This is necessary when geodetic positioning of the receiver 
(cither static or dynamic) is desired, through post-processing the SV data. [igure 3.1 


portrays the flowchart for processing the TI-4100 data. 


A. PRE-PROCESSING 
1. Manipulation of Cassette Recorded Data 
The first step in processing GPS data recorded on cassettes, is to dump the 
raw data from the cassettes to either a 9-track magnetic tape or disk. The Applied 
Research Laboratories (ARL) at Austin, Texas have developed software to accomplish 
this fundamental task. The software at ARL can dump and process the raw data 
recorded on cassettes, and format the GPS data onto a 9-track magnetic tape in the 
SEF (Chapter I].D.). This tape can then be transferred to disk for convenience in 
further processing with a mainframe computer. 
2. Alignment of Common Time-Tagged Data 
The alignment of common time-tagged data items 1s accomplished through use 
of the CON9TR program. 
a. Program Description 
Program CON9OTR, the first of three programs acquired from NSWC 
(Chapter I.), reads the SEF data and formats it into another NSWC specified file. This 
file then becomes the input for the two remaining programs, SOLTOM and KALMN2 
and has all common time-tagged data on the same record. The only calculations 
performed at this stage are variances on the pseudorange data and the delta range 
(Doppler) data which are measures of the signal to noise ratio. 
The program also looks at the “quality vectors” (Texas Instruments, 
1982) of the data on the 9-track SEF tape and returns an integer value representative 
of the data quality at that time line. Four integer values are possible: 
O = good data 
| = good data but this is the first Doppler count 
14 = bad data at that point but there was not a break 
in the Doppler data 
1S = bad data 


11-4100 9-TRACK STANDARD 
CASSETTES EXCHANGE FORMAT TAPE 


CPS DATA FOR GPS DATA 





USER INPUT PROGRAM CONSTR 
FILE 
STATIC | DYNAMIC 
ME SSAGE 
FILE 


USER lhPuT PROGRAH PROGRAH 
FILE 
(SOLUTION TYPE ) 


USER ITMPUT 
FILE 


SOL TOM KAL HA2 
(SOLUTION TYPE) 










RES TOUAL GD OP POSITION 
FILE HIS TOMY 
FILE FILE 


(OPT TONAL ) 


RECEIVER POSITION 
BY $01. UT ION 


RECEIVER POSTT ION 
BY SOLUTION 


TYPE AND TIE TYPE AND TIME 





Figure 3.1) Flowchart for Processing T1-4100 GEOSTAR GPS data. 
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The CON9TR program requires a user input file (Figure 3.1) in addition 
to the GPS data in the Standardized Exchange Format during program execution. 
CONSTR consists of the main program and 22 subroutines. Appendix A gives a brief 
description of each subroutine’s function and nesting. 

b. Input/Output 
The user input file (Figure 3.1) for the CON9TR program consists of the 
following: 
(1) The number of trackers, usually four but may be less than four 
(2) GPS time of the first data point to be put on the new file 
(3) Time to end data on the new file, if the end of the file 
is found before the time to end data, then data will end 
at the end of file 
(4) Range variance bias factor [km2] 
(5) Doppler variance bias factor [km] 
(6) X-coordinate of the initial receiver position that 
will be put on the new files 
(7) Y-coordinate of the initial receiver position that 
will be put on the new files 
(8) Z-coordinate of the initial receiver postion that 


will br put on the new files 


Appendix B gives the user input guide to CON9TR prograin variables and format. 

The first output file from CONSTR (Figure 3.1) writes all common time 
satellite tracking data into the same record block, and Appendix C gives the output 
format. The second output file writes out messages to aid the user in data quality 
analysis, flagging one of four quality vector conditions that exist in the data. 

c. Compution of Variances 
The variances for the pseudorange and Doppler data in the CON9OTR 


program are computed by Equations (3.1) and (3.2), respectively (Texas Instruments, 
1982). 


& t 7 . | 


(0:5 4 Bpyy On 2) ie kopemele 


where the asterick ,*, is used to signify multiplication 


By = code loop noise bandwidth Elz 
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Bpp = code loop predetection bandwidth Hz 
C: = L; carrier power to noise spectral density ratio 
dB-Hz 
N = tracking mode, 1, 2, 3, 4 based on SV mode 
i (chips)? to [m7] conversion: 
(29.305)? for P-code 
(293.05)? for C/A-code 
RK» = variance bias factor [m7] 
@7 51 = [(B, * NVC! * C1 + [0.5* (3.2) 
Bory (101) * K+ Ky 
where, 
PLL = carrier phase-locked loop 
By; = carrier loop bandwidth Hz 
Bop = carrier predetection bandwidth Hz 
= (rad?] to [m-] conversion 
= (c 2nf)* 
c = speed of light 
f = L band carrier frequency 
Wye =0¥.57524 xglO? Hz 
i = 1.2276.x 10? Hz 


; 7 ” 
= variance bias factor [m7] 


7 
£ 
| 


d. Program Modifications at NPS 
The original code from NSWC was written in FORTRAN IV. This 
code was modified at NPS to be compatible with the FORTRAN 77 compilers 
(FORTVS and WF77) on the IBMI-3033 mainframe computer. Changes to the original 
code included the following, 


e Changes in the formatted write statements 


e Changes in variable alignments in the common statements 
(i.e. making sure double precision variables are aligned 


ahead of single precision variables) 
¢ Changes in dimension size of the data arrays 


e Changes in defining variables 


pi 


¢ Changes in the placement of data intialization 


Statements 


e Addition of the Block Data subroutine to contain the 
moved data intialization statements and accompanying 


common statements. 


Executive files, written at NPS to handle the input and output files for 
the CONSTR program, included the auto-double precision commands for all real 
variables using the FORTVS compiler. This modification automatically doubles the 
declared precision in the code and was required due to the use of 60 bit words by the 
CDC Cyber-865 computer and the 32 bit words by the IBM-3033 mainframe computer 
at NPS. 


B. DETERMINATION OF PSEUDORANGE POSITIONS OPO lier 
User receiver position solutions whether fixed (static) or moving (dynamic) may 
be calculated by post-processing the TI-4100 GEOSTAR GPS data. 
l. Static Positioning 
The second program, SOLTOM from NSWC (Chapter I.) computes a position 
solution for a static receiver by using an iterative batch least squares calculation. 
a. Program Description 
Program SOLTOM computes receiver position and time, solving for 
initial corrections to estimates of parameters X, Y, Z, time biases and scaling factor for 
a tropospheric correction. Types of solutions which can be selected through the user 
input file are described in Chapter IIJ.B.1l.e... The program consists of the main 
program and 52 subroutines. Appendix D gives a brief description of each subroutine’s 
function and nesting. 
b. Math Models 
The math models for the batch least squares solution in program 
SOLTOM solve for corrections to the initial estimates of the parameters desired. The 


observation equation is given by Equation (3.3): 


Ov 24V J PN 72 
at lee Xp) ares Y)) a Ae Z5) Wei otG BG DG elie 3 25)) 
where, 
RS = jth range observed from SV of receiver 
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Ao = earth-fixed geocentric X-coordinate of SV 

Yoy = earth-fixed geocentric Y-coordinate of SV 
sy — earth-fixed geocentric Z-coordinate of SV 
Xp = earth-fixed geocentric X-coordinate of receiver 
Yp = earth-fixed geocentric Y-coordinate of receiver 
Ly = earth-fixed geocentric Z-coordinate of receiver 


CB = Clock bias term, SV or receiver depending on user 


input selection 


CD = Clock drift term, SV or receiver depending on user 
input selection 
CA = Clock aging term (frequency drift), SV or receiver 


depending on user input selection 
TR = Tropospheric refraction correction. 


The least squares solution of Equation (3.3) in matrix form, can be written as: 


Bear E = 0 (3.4) 
where, 
B = Matrix containing the partial derivatives of the observation 


equation with respect to the parameters being solved for 


P = Vector containing the differences of the parameter adjusted 
(Pa) and the parameter estimate (Pe), (Pa - Pe) 
E = Vector containing the partial derivatives of the observation 


equations multiplied by the weight matrix, multiplied by the 


data residuals 


The matrix B is obtained from Equation (3.5), 
B= Alwa (3.5) 


je) 


where, 


A = Partial derivatives of data equation with respect to a 
parameter ] 


Al = Transpose of the A-matrix 


w = Weight matrix with can on the diagonal, zeros 


off diagonal, where k = 1| to number of data points, this 


is the variance of noise Ievel on a “k” range measurement 


or Writing Equation 3.5 in summation form, 


Bis = pa OR OP, OR Ty Onnaree (3.6) 


where, 


n = number of data points used in the least squares solution 


Rio = kD gata equation 
ith 


"U 
I 


parameter being solved for 


; jth parameter being solved for 


gc) 
I 


0.2 = ,th variance value 


The A-matrix or design matrix is defined by Equation (3.7): 


Aya OR®: | oP. (ay) 
where, 

R®: = i gata point (3.3) 

P; = jth parameter being solved for 


The inverse of B is the variance-covariance matrix, COV, with the variances 6 of the 


parameter corrections AP that are being solved for, on the diagonal given by Equation 
(3a) 
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5) 


tod 


B! = cov 


there, COV is 2 6 x 6 mairix if six parameters are being solved, or 7 x 7 if the optional 
7th parameter (CA) is ito te determimed. The off diagonal elements are the 
cross-correlations of the paired parameters for that ume line. The vector E in 


summation form is shown in Equation (3.9), 


- 
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E.= Y",_, (éRp°@P.)* 6.7 *R ( 


where, 


Ryo = KD Gara equation 


n = number of data points used in least squares 


calculation 


1 
2 t P : 
~~ = peas weighied variance value 


E = RL,.wA (3.10) 


AP = BIE (3.11) 


and 
AP = corrections io the initial estimaies of unknown parameiers being solved for (X. 


Y. Z, clock biases, and tropospheric refraction correction}. 


The following Equations (3.12) to (3.35) describe ithe various corrections, Time Epoch. 
clock, Ionospheric refraction. Tropospheric refraction. earth rciation. and relativisiic 
effects which are required to correct pseudorange and delta range (Doppler) data 


(Meverhoff, 19$5). 


Time Epoch Correction 


T.., =a (3.12) 


SV C 


where, 


T.. = ume from epoch 
tt = the GPS time-tag for range and delta range (Doppler) 
data 


Je = time of reference for clock corrections 


To correct the T,,, for the end of the week cross over, correct Wey as follows: 


IfT.., < 302400 [s] then, T., = T.y + 604800 [5] 
If T.,, > 302400 [s] then, T,, = T,y - 604800 {s] 


Clock Correction 


Cott) = [Ag + (Ay * Tey) + (Ap * Tey * Tey)l * (3.13) 


AC.\(tt) = [A, + (Ay * {2 * 0) - td) * eee (3.14) 


sv) 
where, 
C.j(tt) = clock correction for pseudorange data 


at GPS time-tag in [km] 


AC,,(tt) = clock correction for delta range 
(Doppler) data at GPS time-tag in [km] 


> 
I 


clock correction bias term [s] 


O 

A, = clock correction drift term [s/s] 

A», = clock correction aging term 
[s/(s)7) 

Te > time fromenoch 
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c = speed of light 299792.458 [km/s] 


td = the Doppler count interval for delta range 


(Doppler) data 


lonospheric Refraction Correction 


Ci (tt) = CP, (tt) - CPo(tt)/(Q,/Q,)* - 1.0 (3.15) 


DC; ,(tt) = (DOP, ,(tt)/Q,) - (DOP, 4(tt)/Q>)] *c (3.16) 
[(Q)/Q>) 7 1.0] Mic 
where, 
C; (tt) = pseudorange ionospheric refraction 


delay in [km] 


tt = GPS time-tag for pseudorange and delta range 


(Doppler) 
CP, = LI pseudorange in [km] 
CP, = L2 pseudorange in [km] 


DOP, , = LI Doppler count 


DOP; 4 = L2 Doppler count 


DC, = delta range (Doppler) frequency shift due 


to ionosphere 


Q, = LI frequency multiplier (154.0) 
Q>5 = L2 frequency multiplier (120.0) 
Vos = Phe nominal satellite oscillator frequency, 10.23 M11z 


In Equation (3.15) the correction, C;,(tt), is used to correct the pseudorange that ts 
time-tagged, tt. In Equation (3.16) the correction, DC;,(tt), is applied to the delta 


range (Doppler) data to correct it for time-tag, tt. 


p47) 


Tropospheric Refraction Correction 
The Tropospheric refraction computation can be divided into two parts, 
weather dependent and receiver geometry calculations (Meyerhoff, 1985). 
Weather Dependent Correction 


First, convert temperature to degrees Kelvin and convert humidity to a 


fraction: 
ie. 260 eal) 
RH = H,/100.0 (3.18) 
where, 
T, = air temperature in degrees centigrade 


Hy = the relative humidity in percent. 
Second, compute the surface partial pressure, E,, from temperature and relative 
humidity. 


E, = RH * 35.65 * 10 exp[7.617-(2285.0/T,)} (3.19) 


Third, the components of the zenith value of tropospheric delay, Zary and Zo, are 


calculated from, E,,, and atmospheric pressure. 


Zary = 2.276 * Pb * 0.01 (3.20) 
Zwet = (470 * E,exp(1.23)/T,) + (1.705 * 10° (3.21) 
* A, * E,exp(1.46/T;,) 
where, 


exp = exponent 
Pb = the atmospheric pressure in [mb] 
Ay = temperature lapse rate set at 0.006 °C/m 
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exe — Exponent it should be noted that Lary and Z,.,., need only be recomputed, 


when the weather data changes. 


Receiver Geometry Dependent Correction 
First, to find the elevation angle, EV(t), of the satellite at time, t, take 
the inverse sine of the dot product of the unit vector from the receiver to the satellite 


and the unit position vector of the receiver. 


a 
R 





ee: Ree 
E(t) — sin RS (3522) 
Gr vie 
Mere, 
aia 
Ie Vector fromthe receiver to the satellite 
MG,, = vector’s magnitude 
—_—_ eo @ * 
R = position vector of the receiver 
——— 
MG, = magnitude of R 
MG... = [(X(t)- X,)* + (Y(t) - Y,)* + (Z(t) - Z,)*4!/4 628) 
= 2 2 ali? | 
ice [x.> eyeerez | (3.24) 
where, 
Se Bes Ee = the estimate of the receiver's 
position in an earth-fixed 
geocentric coordinate system 
rewriting Equation (3.23), EV(t) can be written as, 
ev — 1Cx(t)- X,)* A,) (3.25) 


sme) ea) ZOMG... * MG,) 
Alene: 
X(t), Y(t), Z(t) = the position coordinates in a earth- 


oy, 


fixed geocentric coordinate system 


Or 
EV(t) = sin’! [SEV(t)] (3.26) 


Second, the multiplication factors for dry and wet as a function of elevation are 


calculated from the elevation angle, EV,at GPS time-tag (tt). 


Fgry(t) = 1/ [sin EV(t)] + {0.00143/[(tan EV(t)) + 0.0445} (3.27) 


Fye(t) = 1/ [sin EV(t)] + {0.00035/[(tan EV(t))+ 0.17} (3.28) 


Third dry(t - td) and F\,,(t - td) are calculated in a similar manner as Equations 
(3.27 and 3.28). 


Total Tropospheric Refraction Correction 


TR(tt) = (Zary * Fg) peer (3.29) 


DTR(tt) = (Zgry * [Fury = Baye td) ee (3.30) 


F wet Es td)}} 
where, 


TR (tt) and DTR(tt) are range corrections are in [km]. 


Earth Rotation Correction 
An earth rotation correction is required due to the rotation of the earth 


during the signal propagation. 


Co, — WYO* XCD Ve (3.31) 
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DC, = [({¥(t) - Y(t - td)} * X,) - ((X(t)-X(t - td)} * YJ] * RRE/c (3.32) 


where, 

X(t), Y(t), Z(t) = position coordinates in an earth- 
fixed geocentric cartesian 
coordinate system at pseudorange 
transmit time “t” used with data 


time tagged “tt” 


GP — = estimate of the X and Y coordinate 
in an earth-fixed geocentric 
cartesian coordinate system 
RRE = Rotational rate of the Earth corrected for motion 
of the equinox 


(0.00007292115855 [rads/s] 


td = Doppler count interval for delta range 


(Doppler) data. 
Relativistic Correction 
The correction for the relativistic effects is due to the eccentricity of the 

satellite’s orbit. 

Rel(tt) = (2 * {GM]!/2/c) *e * A * sin [Ek(t)] (3.33) 
at the beginning of the Doppler count, 

DRel(tt) = (2 * {GM]!/2/c) *e * A * sin [Ek(t - td)} (3.34) 
which yields, 

DRC(tt) = RC(t) - RC(t -td) (25) 
where, 
tt = GPS time-tag for both pseudorange and delta range 


(Doppler) data 


3] 


td = Doppler count interval for delta range (Doppler) 
data 


t = GPS time-tag at transmission time for data time- 


tagged tt 
e = eccentricity of the satellite orbit 


A = the square root of the semi-major axis of the 


satellite’s orbit 


Ek(t) = eccentric anomaly for the satellite at pseudorange 


transmission time t 


Ek(t - td)= eccentric anomaly for the satellite at the 


beginning of the Doppler count interval 


GM = WGS 72 value of the Earth’s gravitational constant including 
atmosphere 398600.8 km>/s* 


RC = computed range 


these are the defined constants and variables. 
c. Input/Output 

There are two input files for the Soltom program (Figure 3.1). The first 
file is the user input file, where the user sets which types of solutions will be calculated 
and which data solutions will be printed out. The user input guide for this file is 
described in Appendix E. The second input file for SOLTOM is an output file from 
the CONSTR program containing records that are aligned with common time-tagged 
GPS data items. 

There are also two output files for SOLTOM (Figure 3.1). The first file 
is the updated receiver position, time and other print options such as solution sigmas 
and clock bias terms, (as determined by the user input file). The second output file is a 
data corrections residual file, which can be used in plotting routines and for input to 
other software that handle multiple station positioning (Meyerholf, 1986). 

d. Program Modifications at NPS 
The source code from NSWC for the SOLTOM program was written in FORTRAN 


IV which has been modified to run at NPS on the IBM-3033 mainframe computer. 
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Similar changes to the source code in the CON9TR program were also made to the 
SOLTOM program including the executive files (Chapter III.A.2.).. The original code 
was Written to be able to handle the precise ephemeris data as well as the onboard 
broadcast ephemeris data in the navigation data message from the SV. The test data 
for SOLTOM was generated at NSWC using the Broadcast ephemeris only, and as 
such, the modified program version at NPS has the Precise Ephemeris subroutines and 
corresponding, call, statements commented out. 
e. Solution Techniques Through Least Squares 

There are ten types of solutions possible with the SOLTOM program 
which are dependent on the user input file (Meyerhoff, 1985). Appendix E gives the 
details for the tvpe of solution desired. Of these ten types, nine are batch least squares 
solutions and one solution type five is a Newton-Raphson solution (a geometrical 
solution, not a least squares fit). The type-five solution can be used for either static or 
dynamic positioning performing an independent solution at each time line. The nine 
batch least squares solution types fall into three data classes, Range (pseudorange), 
Range Difference (delta range Doppler), and Bias Range. 

The reader is IVARNED, that the Bias Range (solution tvpe ten) has still 
not been tested by NSWC, and thus it should not be used at NPS. Solution types one, 
four, six, and nine are Range (pseudorange) solutions. Solution types two and seven 
are Range Difference (Doppler) solutions. Finally, solution types three and eight are 
combination solutions of the Range (pseudorange) and Range Difference (Doppler) 
data. 

The SOLTOM program has the capability to run more than one 
iteration of batch least squares, if needed (e.g. for bad initial receiver position or time 
bias problems). The user has control via the user input file to set how often the least 
Squares solution will be calculated and printed out. When a solution is to be 
computed, the program takes the B-matrix containing the partial derivatives of data 
equations with respect to the parameters being solved for, and multiplies by the 
variances of the range measurements (Equations 3.5 and 3.6), inverts the B-matrix, 
which is the variance-covariance matrix with variances of the corrections AP on the 
diagonal, and computes the corrections to the initial estimates of the parameters being 
sought, AP (Equation 3.11). Then Equation (3.36) gives the updated station position 


coordinates, clock bias, drift, aging, and Tropospheric refraction correction: 


a ear. (3.36) 
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Y= ¥ + APY 
Z=Z+ AP, 
A SN 


Ajt = Ajt mE ai 
Ds? 
a EY fo 


where, 
X, Y, Z = updated earth-fixed geocentric cartesian 


coordinates through batch least squares 


Ag = clock correction bias term in [s] 


A, = clock correction drift term in [s/s] 


A» = clock correction aging term in 
[s/(s)7] 


(optional 7th parameter to solve for) 
TR = tropospheric refraction correction 


The new solution takes all previous solutions and averages up to the current time line, 
before producing the iterated solution. 
2. Dynamic Positioning 
The third program, KALMN?2 from NSWC computes a position solution for a 
moving (dynamic) receiver, using the UDU’ factorization of an eight-state extended 
Kalman filter (Chapter [.). 
a. Program Description 
Program KALMN2 uses the first four-states for receiver position (X, Y, 
Z) and receiver clock bias, and the other four-states for receiver velocity along three 
axes (X, Y, Z) and receiver clock drift. By using the Doppler phase (change in range) 
measurements to smooth the pseudoranges, signal multipath and measurement noise 


effects can be suppressed significantly to produce a viable solution (Meyerholff and 
Evans, 1986). 
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Note: The four-states of the velocity and clock drift arc still under 
testing at NSWC. 

The KALMN2 program consists of the main program and 43 
subroutines. Appendix F, gives a brief description of each subroutine’s function and 
nesting. 

b. Math Models 

The math model for the UDU’ factorization of the inverse of the 
B-matrix (which contains the partial derivatives of the observation equations with 
respect to the parameters being solved for) used in program KALMN2 is described in 
Bierman (1977) and Meyerhoff (1986). 

The TI-4100 receiver provides two pseudorange and two Doppler 
measurements (via LI and L2 frequencies), for as many as four satellites 
simultaneously. The following Equations (3.37 to 3.42) describe the various corrcctions 
and conversions; receiver clock bias, continuous count (Doppler) to interval (Doppler) 
data, and lonospheric refraction, which are required for the Doppler and the 


pseudorange data (for the Ionospheric correction) (Meyerhoff, 1985): 


Receiver Clock Bias Correction 
tt = ttr- RCB (3738) 


Winere, 
ttr = receiver clocks GPS time-tag for range and delta 


range (Doppler) data 
tt = GPS time-tag correctcd for receiver clock bias 


RCB = cstimate of the receiver clock bias which 1s 
initially set to zero and an improved estimatc 


of “RCB” is made aftcr each time line 


The time delay data are multiplied by the speed of hght “c” and 


converted to a pseudorange in [km]. 


Continuous Count to Interval (Doppler) Data Conversion 
DOr tee = DC) (tee DCHL - tay OSE 1) (3.38) 


ae 


DOP, 5(tt); = DC4(tt), - DCa(tt - td), - (td * OF FSET2) (3) 
where, 
td = Doppler count interval for delta range (Doppler) data 


DC (tt - td), DCy(tt - td)= Li and L2 continuous 


Doppler count data at the delta range interval 


DC, (tt) , DC,(tt) = Ll and L2 continuous count Doppler 
data at the end of the delta range interval 

OFFSET, , OFFSET, = Doppler offsets used by the receiver 
for LI and L2, -6000 and +7600 respectively 


DOP, ,(tt), DOP, (tt) = LI and L2 interval Doppler data 


for time-tag “tt” 
1 = satellite from which data originates 


The delta range (Doppler) data is obtained by converting the LI interval Doppler data 


for each satellite. 
RD(tt); - (-DOrT (tu y =. ¢ Q,) ae (3.40) 


where, 


V 


os = nominal satellite oscillator frequency, 10.23 MHz 


Q, = LI frequency multiplier (154) 
RD(tt) = delta range (Doppler) for one satellite in [km] 


lonospheric Refraction Correction 


RCI(tt), = CPI(tt), + C1O(tt); (3.41) 
RDCl(tt); = RD(tt); - DCIO(tt); (3.42) 
where, 
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CPI = uncorrected pseudorange measurement in [km] 


RD(tt) = uncorrected delta (Doppler) range measurement 


ClO(tt) = tonospheric refraction correction in [km] 


for pseudorange data (refer to Equation 3.15) 


DCIO(tt) = ionospheric refraction corrections in [km] 
for delta (Doppler) range data (refer to Equation 
Sane 


Smoothing For Noise 

The pseudorange data are contaminated due to signal multipath and 
other noise sources (Meyerhoff and Evans, 1986). To suppress the noise duc to 
multipath in the pseudorange data, they are smoothed by the delta range (Doppler) 
data. The next set of Equations (3.43, 3.44, and 3.45) describes this smoothing: 


BN Gt) = Kitt) - DR(tt,ttfst) (3.43) 
papel DBi(tt)/N | (3.44) 
eu) = Baygltt) + DR(tt,ttfst) (3.45) 
where, 


tt = GPS time-tag for pseudorange and delta range (Doppler) data 


ttfst = GPS time-tag for the first good pseudorange and 


Doppler count of this data span 
R(tt) = Ionospheric corrected pseudorange in [km] 
DR(tt,ttfst) = lonospheric corrected delta range (Doppler) data 


N = number of data points in the span 


NX Nth estimate of pscudorange at time “ttfst” 


a0 


Baye = average value of the first pseudorange 


Ro w(tt) = smooth pseudorange, if the absolute value of 
the difference between the raw pseudorange and 
smoothed pseudorange exceeds 20 [m], the 


data point is considered bad. 


Both the pseudorange and delta range (Doppler) data must be corrected for Time 
Epoch, clock errors, earth rotation, relativistic effects, and Tropospheric refraction. 
These formulae are the same ones used for the iterative least squares solution in the 
SOLTOM program (Equations 3.12, 3.13, 3:14, 3:31, 3.32,°3°33)323ameo53 eee ee 


The corrected pseudorange 1s given by, 
RC(tt), = Ro (tt): + Coj(tt); - Rel(it). - [TR(tt): Au] gee GRO eaten oe (3.46) 


RDC(tt) is calculated in a similar manner as (3.46) for the delta range (Doppler) data. 


The following Equations (3.47 to 3.51) describe computed ranges and residuals for 


pseudorange and delta range (Doppler) data (Meyerhoff, 1985): 


Computed Ranges 

CR(tt); = [(X(); - KX)? + (¥(O-Yr)* + (Z(0; - Z4}4/4 (3.47) 
for the range at the beginning of the Doppler interval, 

CER (tt 10d) eo 2) as [tea td)}}2 tA (U2) are td)}}7 (3.48) 


+ {Z(t2); - [Zr - (ZVr * td} 2p? 
where, 
X(t), Y(t), Z(t) = position of i“) satellite in an earth- 
fixed geocentric cartesian coordinate 


system at transit time t 


Xr, Yr, Zr = estimate of the receiver’s position in an 
earth-fixed geocentric cartesian coordinate 


system 
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X(t2), Y(t2), Z(t2) = position of it} satellite in an 


earth-fixed geocentric 


cartesian coordinate system at 
time t2 


XVr, YVr, ZVr = estimate of receiver’s velocity 


tt = GPS time-tag corrected for receiver clock bias 


td 


Doppler count interval for delta range (Doppler) 
data 


Therefore the computed range difference for the individual satellites 1s, 


CRD(tt): = CR(tt); - CR(tt - td); (3.49) 
Residuals 


The pseudorange residuals are given by, 


pega (tt) = sett) C(t) REB 


(3.50) 
and for the delta range (Doppler) data, 


RDomc(tt); = RDC(tt); - CRD(tt), - RCD (3.51) 


where, 


KS CHIGOE 


corrected pseudorange (Equation 3.46) 


RDC(tt); = corrected delta range (Doppler) 


CR(tt); = computed pseudorange 


DCR(tt); = computed delta range (Doppler) 


RCB estimate of the receiver clock bias 


RCD 


estimate of frequency bias of the receiver oscillator 


which is also receiver clock drift. 


a9 


tt = GPS time-tag corrected for receiver clock bias 


td = Doppler count interval for delta range (Doppler) 
data 


ith 


— 


satellite being tracked 


A state vector SP is corrected and improved over time, 1s composed of 
the receiver position in an earth-fixed geocentric cartesian coordinate system (Xr, Yr, 
Zr) and receiver clock bias (RCB). The last set of Equations (3.52 to 3.58) describes 
the observation equation and development of the state vector for the pseudorange data 
(Meyerhoff, 1985): 


Observation Equation 

D(tt)) = CR) aan eb (3.52) 
or rewriting it as, 

D(tt); = [(X(t); - Xn)? + (Y(v, - Yr)? + (Z(t); - Z)7]!/4 + RCB (3.53) 


The H-matrix is defined as containing the partial derivatives of the 


observation equation with respect to the state vector: 


Uh es OD(tt)/OXr = -(X(t), - Xr)/CR(tt); (3.54) 
Hj = ODUD/OYr = (x(t) a) Gite 

= oD(tt)/oZr = -(Z(t); - Zr)/CR(tt); 

Hg | = 0D(tt//ORCB = 1 


The PA-matrix is defined as a 4x4 variance-covariance matrix for 


— 
—_— 
os) 
“a 
p= © 


position and receiver clock bias, with off diagonal elements zeroed out and the diagonal 
elements set to the initial variance of Xr, Yr, Zr, and RCB. The Kalman vector K is 
given by the matrix equation, 


K(tt) = PA(t-) * 1H. *(H, * PA(t-) * HT, + Rey! a2) 


where, 
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PA(tt-) = variance-covariance matrix at time tt before it is 


updated for this data point 


R = variance on the pseudorange observation 


i = iD satellite being tracked 


Once the PA-matrix is updated it is defined by, 


Battet )) SRA) eek (tt) eM; * PA(it-) (959) 


State Vector 


The final two equations define the updated state vector SP(tt+): 
Srait+ ) > SP(tt-) + K(tt); * Rome(tt); (2257) 
or 


i) — Na) Kj (tt). * Romic(tt): (3.58) 
peer p= ya(-) + KA(ct))* Romce(tt); 

ZG oe Ze) Ka (tt) = Romic( tt): 
RCB(+) = RCB(-) + Ka(tt); * Rome(tt); 


where (+) means the updated state vector and (-) means before the updated state 
vector. 
c. Input/Output 

There are two input files for the KALMN2 program (Figure 3.1). The 
first file is an output file from the CON9TR program which contains records that are 
aligned with common time-tagged GPS data items. The second file is a user input file 
Which sets the types of solutions and/or constraints on that solution, along with which 
data solutions will be printed out. The user input guide for this file 1s in (Appendix G). 

Mere are titee possibie output files for KAIANEN2 (ligure 3.1). The 
first file is the updated receiver position and time and any other print options that are 
set by the user input file. The second file is a position history file that 1s used in 
printing out the updated receiver position and time file. The third file is an optional 


GDOP file, that can be set in the user input file. 


4| 


d. Program Modifications at NPS 

The source code from NSWC for the KALMN2 program was written in 
FORTRAN IV which has been modified to run at NPS on the IBM-3033 mainframe 
computer. Similar changes to the source code for the CON9TR and SOLTOM 
program were also made to the KALMN?2 program including the executive files 
(Chapter III.A.2.).. The original code was written to handle both the precise as well as 
the onboard broadcast ephemeris data in the navigation data message from the SV. 
The test data for KALMN2 was generated at NSWC using the broadcast ephemeris 
only, and as such, the modified program version at NPS has the precise ephemeris 
subroutines and corresponding, call, statements commented out. 

e. Solution Techniques Through Kalman Filter 

Using a UDU’ factorization of an eight-state extended Kalman filter 
(with smoothed pseudorange measurements) in program KALMN2, there are nine 
types of solutions that may be calculated depending on the user input file (Appendix 
G). 

Note that, as was stated previously in Chapter III.B.2.a., the velocity 
solutions from the four-velocity states are still being tested at NSWC. 

The initial receiver position (X, Y, Z) and clock biases are updated using 
the corrected state vector from the Kalman filter (Equations 3.57 and 3.58). If for a 
given time line the absolute value for the ionospheric correction CIO(tt) exceeds 20 m, 
that data is considered bad and flagged as such, with a quality vector set to 15. [eight 
constrained (especially in marine areas) allows the user’s receiver to track four satellites 


with poor GDOP and still be able to produce a good solution. 
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IV. TEST RUN RESULTS AND VALIDATION 


The testing of the three programs CONSTR, SOLTOM, KALMN2 was 
accomplished by using two independent data sets in the SEF on 9-track magnetic tape 
as received from NSWC. Included also were hardcopy printouts, listing the solutions 
from those data sets using the CON9TR, SOLTOM, and KALMN2 programs on 
NSWC’s CDC Cyber-865 computer. A direct comparison of the output listing files 
containing the receiver position solutions and times for the SOLTOM and KALMN2 
programs from both the NSWC’s CDC Cyber-865 and the NPS IBM-3033 computers 
could then be made. 

After the original source codes and some algorithms were modified for 
compatibility on the NPS IBM-3033, the programs were test run, using the FORTVS 


compiler with an auto-double precision option discussed previously (Chapter II].A.2.). 


eo were Os| TIONSORUTION RESULTS 

Table II shows the differences between results computed at NPS and NSWC 
using the NSWC test data set (from Henderson Point, 85006) for the SOLTOM 
program. Running the CONSTR and SOLTOM programs, position solution results 
(X, Y, Z) were duplicated to the millimeter level for solution type seven; to the 
centimeter level for solution types one and two; to the decimeter level for solution 
types four and six; and to the meter level for solution type nine (Appendix E). 

The above differences in the SOLTOM solutions are considered insignificant due 
to the hmitations of the observed data (Meyerhoff, 1986). 


Pee Nave POSITION SOLUTION RESULTS 

Table III shows the differences between results (beginning and ending solutions 
using the test data set) computed at NPS and NSWC using the NSWC test data (from 
FHlenderson Point, 95005) for the KALMN2 program. Running the CONSTR and 
KALMN2 programs, position solutions and clock biases were duplicated to the 
millimeter and 10°© second level respectively. 

It should be noted that not all options (Chapter II1.B.2 and Appendix G) within 
the programs have been checked and tested. At this point it can be said that these 
programs have been established and validated on the IBM-3033 at NPS to the degree 


of being able duplicate NSWC’s results if the same user options are followed. 
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TABLE II 
CONSTR AND SOLTOM SOLUTION RESULTS 


Difference Between NPS and NSWC Solutions 


F X,Y,Z . Number | 
50106 10n peace “(286 Bias Clock Drift Of Data Points 
Type (kilometers) (seconds) (seconds) , 


oh -5 -9 -11 
Range Data 10 10 10 
Receiver Clock 


WO -5 -11 
popalee Data 10 10 
Receiver Clock 


FOUR -4 
Common Range 10 
Receiver Clock 


S/1X =i 
Range Data ue 
SV Clock 


<popoler Mee ikem 


NINE = 
Common Range 10 
SM 1 | OE 
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ibe I 
CONSTR AND KALMN2 SOLUTION RESULTS 


Difference Between NPS and NSWC Solutions 
GPS Time-Tag 1 Ve) Clock Bias Number of 
(kilometers) 


(seconds) (seconds) Doppler Edited 
Points 


First Solution 

14129. 00 107° 1079 N/A 
Last Solution 

17499. 00 107° 1078 





45 


V. SUMMARY 


Ilow to process the GPS satellite data from raw cassettes off the TI-4100 
GEOSTAR receiver to a reduced position and time solution has been discussed. 

The capability of processing raw data from cassettes to SEF at NPS 1s a task for 
the atu: 

NPS at this time has the capability, starting with a SEF 9-track magnetic tape, to 
use the CON9TR program (where the common time-tagged data are aligned), to select 
either the SOLTOM (an iterative least squares for positioning a static receiver) or the 
KALMN2? program (which uses a Kalman filter for positioning a dynamic receiver) and 
to compute position and time solutions. The calculations for program CON9YTR, and 
the math models for programs SOLTOM and KALMN2 were documented. The 
changes required to modify the original codes and algorithms for these programs 
written for NSWC’s CDC Cyber-865 computer to make them compatible with the 
IBM-3033 computer at NPS have been described. Testing of the programs shows a 
successful conversion (with the restrictions discussed in Chapter IV.) of the programs 
to the IBM-3033 mainframe computer. This software forms a basis for a GPS 


programs library that will aid research using the GPS at NPS. 
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VI. RECOMMENDATIONS 


Since the conversion of the CONSTR, SOLTOM, and KALMN2 programs has 
been successfully tested using the two independent data sets from NSWC, it 1s 
recommended to explore other options within the programs set by the user input files 
for the SOLTOM and KALMN2 programs. Options that should be explored for the 
SOLTOM program include the use of the precise ephemeris and clock files, the use of 
more than one iteration for the solutions, and the calculation of Newton-Raphson 
range solutions followed by a comparison with the iterative batch least squares 
solutions for a given time line. The KALMN2 program should be explored by using 
the precise ephemeris and clock files, turning the height constraint on and off in the 
user input file and comparing the solution results. The refraction correction should be 
turned on and off to see how the solutions are affected (this is another option in the 
user input file). 

It is further recommended to explore the possiblities of rewriting or adding codes 
and algorithms to improve the solution corrections for updating the four-state vectors 
of the velocity in the extended eight-state Kalman filter making them reliable for use in 
velocity solutions. 

In the future TI-4100 software for carrier phase data, differential positioning etc., 
should be acquired and modified to be compatible with the IBM-3033 computer at 
INPS. 
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APPENDIX A 
CONSTR SUBROUTINE DESCRIPTIONS AND NESTING 


INP y= This routine reads the user inputs. 
GETDAT - This routine reads the next set of satellite 


tracking data from the 9-track Standard tape 
and also updates weather, Navigation data, 
almanac and calibration. 


RDDTL - This routine initializes the quality vector 
lags and controls branching of the program 
depending on a particular crror flag. 


ERSGES- This routine controls branching SS toh irc 
for several error and terminating conditions. 

UDNTDT - This routine USS non-tracking data so that 
it can be used, alter it determines wien we 


of non-tracking data was read from the 9-track 
Standard tape. 


PRES: This routine prints tape solutions. 


UDWITDT - This routine gets the weather data from the 9-track 
Standard tape and stores it in common that will be 
used to calculate the tropospheric correction. 


UDSATD - This routine reads the 9-track Standard tape and 
defines new satellite orbit clock paramcters. 


FDTRN - This routine will give the tracker number for the 
satellite that the new data 1s for. 

DEDPOF - This routine stores the values for the LI and L2 
Doppler offsets so that they may be used in the 
program. 

PRTPDT - This routine prints any unused 9-track Standard 
tape data. 

UDTRD - This routine decides if the tracking data 1s 


continuous count Doppler range or data quality 
measurement blocks and then Calls the correct 
subroutine to handle the correct data type. 


DEQUDT - This routine defines data qualitv and then stores 
the data quality information. The satellite 
tracker, mode and the quality vectors are used to 
determine if the data are good. The remaining data 
values in this block are used to calculate errors 
on the tracking data. 


DERGDT - This routine stores the range data and the signal 
to noise data in correct data array that will be 
processed. 

DEDPDT - This routine converts the accumulated Doppler 


hase data mto range diflerence data. The 
doppler phase data is corrected for offset, and 
then it,is differenced with preceding data if the 
preceding data was good. 
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GNINP - 


GNID - 


INSTR - 


BLDATA - 


SDAT - 


FSIG - 


This routine reads the logical unit number of the 
input file and returns the information read one 
item at a time. 


This routine controls how data items are read by 
reading a four digit data tvpe number. 


This routine initializes the read package of the 
9-track Standard tape and sets up a data block 
for the data items read off of the tape. 


This routine contains the data that has been read 
off of the 9-track Standard tape. 


This routine takes the site name and the data 
interval from the 9-track Standard tape read 
common and passes thesc values. 


This routine determines the difference between 
the raw range and the range diflcrence data for 
the satellite being tracked. 


BLOCK DATA - This routine was added at NPS to make this 


program compatible with the FORTRAN 77 
conipilers on_the [BM-3033 mainframe computer. 
It consists of common statements and data 
statements that the rigorous compilers 

at NPS would not permit where thev were placed 
in the original source code of CORDTR. 


MAIN PROGRAM CONSTR 
INPUT 
GETDAT 
RDDTL 
ERITR 
UDNTDT 
PRTSL 
UDWTDT 
UDSATD 
~_ FDTRN 
DEDPOF 
. PRIPDT 
UDTRD 
cu — 
FDTRN 
DERGbL 
FDTRN 
DESP 
. FDTRN 
GNINP 
. GNID 
INSTR 
BLDATA 
GNINP 
RDDTL 
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APPENDIX B 
INPUT GUIDE FOR PROGRAM CONSTR 


( by S. L. Meyerhoff, NSWC) 


This is the user_input file for the program with variable names and descriptions. The 
input is format free. 


line | 
1: NSAT = number of trackers (4) 
2; ISS = GPS time of first data point to put on 
new file 
Bf LES = time to end data on new file, if EOF 
found before TES then data will end at 
Bor 
4: K2 = range variance bias factor (kilometer) 
5: K4 = doppler variance bias factor (kilometer) 
line 2 
ioe = X coordinate of initial receiver position 
that will be put on the new files 
PEG = Y coordinate of initial receiver position 
that will be put on the new files 
Suez. = Z, coordinate of initial receiver position 


that will be put on the new files 


SO 


APPENDIX C 
TEXAS INSTRUMENTS GPS SATELLITE RECEIVER FILE FORMAT 


( by S. L. Meyerhoff, 1984, NSWC ) 
Data from the DMA version of the TI-4100 GPS satellite receiver can be put 
into the following format. In this format all common time satellite tracking data will 


be in the same record block. The format for this file is as follows. 
var type meaning 
Record One 


word | ITYPE intg This is the data type for 


the next record 
Record Iwo 
If [TYPE equals one then record two will contain initial information. 


word § “RCTOrr real frequency bias offset for LL] 


word 2 RL2ZOFF real frequency bias offset for L2 
word 3 PB real pressure in (millibars) 

word 4 TP real temperature (degrees Celsius) 
word 5 HD real relative humidity (percent) 
word6 X real X coordinate of initial 


receiver position (kilometers) 
word 7 Y real Y coordinate of initial 

receiver position (kilometers) 
word 8 Z real Z coordinate of initial 

receiver position (kilometers) 
WiOTG. 75 ed in intg maximum nuimber of trackers used 


will be 4 or larger 


word 10 DIREC real user given data interval 
(seconds) 

word I1 SITE(1) char site identification for data 

end eos! Vit 2) char site identification for data 


Sl 


(continued) 


If ITYPE equals two then record two will contain satellite ephemeris data (navigation 


data message) 


word 1 J intg tracker number 

word 2 ISAT(J) intg satellite prn number 

word 3 AS(1) real satellite clock bias in seconds 

word 4 AS(2) real satellite clock drift 
(seconds/second) 

word 5 AS(3) real satellite clock aging 
(seconds/second*second) 

word 6 CICS real cross track cosine amplitude 
(radians) 

word 7 CISS real cross track sine amplitude 
(radians) 

word 8 CRCS real radial cosine amplitude 
(kilometers) 

word 9 CRSS real radial sine amplitude 
(kilometers) 

word 10 CUCS real in track cosine amplitude 
(radians) 

word 11 CUSS real in track sine amplitude 
(radians) 

word 12 DNS real corrections to mean niotion 
(radians) 

word 13 ES real eccentricity 

word 14 IDOTS real inclination rate 
(radians/second) 

word 15 IOS real iiclination (radians) 

word 16 MOS real mean anomaly (radians) 

word 17 OMEDS real night ascension rate 
(radians/second) 

word 18 OMEGS real right ascension (radians) 

word 19 SQAS real square root of semi-major 


axis (kilometers) 
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word 20 ADE real age of data ephemeris 


(seconds) 
word 21 TOCS real clock epoch (seconds) 
word 22 TOES real ephemeris epoch (seconds) 
word 23 IWKNOS intg week number 
word 24 WS real argument of perigee 
(radians) 
word 25 IEDATS intg sv health 
word 26 ADC real age of data clock (seconds) 


If TYPE equals three then record two will contain the tracking information. 


word 1 ay real GPS time-tag for data 
that follows (seconds) 
word 2 ISAT(1) MiesekeNeruinber for the 
to to intg satellite being tracked 
word NIR+1 ISAT(NTR) On trackers 1 toOeNTR 
wordsyiR+2  CRIVT) The L1 pseudorange on 
to to real each tracker (seconds) 


word NTR*2+1 CRI(NTR) 


word 1 R*2+2 CR2(1) The L2 pseudorange on 
to to real each tracker 

monde k*3+1 CRZCUNTR) (seconds) 
WordNWTR 3+2 DOP!) The LI continuous 

to to real count doppler data on 
worden h*4+ 1 DOPRTGNTI<) each tracker (cycles) 
WOr@NIR*4+2 D@P2(1) The L2 continuous 

to to real count doppler data on 
worcmes Rh’ Ot! DOPZENTR) each tracker (cycles) 
wonder R*5+2 SGRI(1) iiiie variance on the L1 

to to real pseudorange 
WOordMNTR*6+1 SGRI@NTR) (kilometer*kilometer) 
werd NTR*6+2 SGR2(1) ive variance on the LZ 
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to to real pscudorange 


word NTR*7+1 SGR2(NTR) (kilometer*kilometer) 

word NTR*7+2 SGDI(NTR) The variance on the L1 
to to real doppler 

word NTR*8+1 SGDI(NTR) (kilometer*kilometer) 

word N#R*38 +2 SGDzZap The variance on the L2 
to to real doppler 

word NTR*9+1 SGD2(NTR) (kilometer*kilometer) 

word NTR*9+2 MQVEL(1) The quality vector for 
to to intg all trackers 

word NTR*10+ 1 MQVEL(NTR) 1S = loss of lock 


14 = bad data, no 
loss of lock 

1 = good range data 
bad doppler data 

O = all data good 
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APPENDIX D 


SOLTOM SUBROUTINE DESCRIPTIONS AND NESTING 


GAIN DIO - 


PEXYZ - 


SIGAYZ - 


MATRIX - 


LINV2F - 


VaIViL EP - 


VMULFM - 


INDAT - 


NEWSOL - 


LSTSOL - 


GETDAT - 


Pewrpl - 


DEFBR - 
CRDAT - 


BSADR - 


CISDAT 2 = 


ION - 


LORTK ID - 


This routine reads the ee for program SOLTOM 
iroM tie Output OM_CONITEK, 


This routine converts position in the geodetic. 
coordinate system to the earth-fixed geocentric 
coordinate system. 


This routine converts sigmas in the Geter 
system into sigmas in the earth-fixe 
geocentric coordinate system. 


This routine performs matrix operations by 
calling routines in the IMSL library. 


This routine is in the IMSL and does a matrix 
inversion. 


This routine is in the IMSL and does matrix 
multiplication. 


This routine.is in the IMSL and does a 
transpose of matrix multiphcation. 


This routine will pass initial values to some 
data arrays. 


This routine sets up the configuration for a 
for a given solution routine. 


This routine calculates a sequential least 
pies solution for the receiver position and 
CIO. 


This routine reads the next set of satellite . 
tracking data from CON9TR output file and it 
will also update weather, Navigation data, 
almanac and calibration. 


This routine converts the accumulated Doppler 
phase data into range difference data by 
correcting phase for offset and then 

differencing it with data from before 

(if before data 1s good). 


This routine defines bias range 


his routine delines correction lor receiver 
1as. 

This routine will correct data for receiver 

clock bias, drift and also correct GPS time 

tag for the end of week cross over. 


This routine converts data into range 
residuals and corrects data [or crror sources. 


This routine corrects for ionospheric 
correction. 


This routine corrects range difference data. 
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EPHDAT - This routine defines the following values 
from the ephemeris data, satellite tracker 
number, satellite bias, drift and aging terms, 
semi-major axis of satellite orbit, 
eccentricity of the satellite orbit and epoch 
time for clock terms. 


CLOCK - This routine adds satellite and clock error 
and converts it to range measurements In 
kilometers and also performs the same for range 
difference data. 


SATPOS - This routine will find the position of the 
satellite from the Navigation Message or from 
the Precise Ephemens: 


ESTPOS - This routine finds the position of the 
satellite from the Precise Ephemeris data. 
ZANG - This routine finds the Zenith angle by doting 


the \aee ios ious of the satellite receiver 
osition vector to the vector from the receiver 
o the satellite. 


RES” - This routine finds the residual by calculating 
the range to a satellite and subtracting it 
from the observed range. This is also 
performed for range difference data. 


SMOT - This routine determines an earth rotation 
correction. 

TROP - This routine finds tropospheric refraction 
correction. 

REL - This routine finds a relativistic correction 
to the data. 

FSIG - This routine finds the variance (sigma squared) 


for the raw range and range difference data for 
the satellite being tracked. 


SOLTNS - This routine calls the different routines to. 
sum data in the correct part of the “B” matrix 
and “E” vector depending on which options have 
been set by the user for type of solution. 


SMPTI - This routine calls a routine to define additions 
to “B° matnix and E vector pilic direc een 
called routine on where to store this 
mformation for the “B” matrix and “E” vector 
(for range data only). 


UDRRC - ea range solution matrix receiver clocks. 
Then store data i matrix so that solution can 
coinputed. 

SMPT2 - This routine calls a routine to define 


additions to “B” matrix and “E” vector plus 
directing the called routine on where to store 
this information for “B” matrix and “E” vector 
(for range difference data only) 


UDRDRGe Update range difference solution matrix 
recciver bias. Thus routine stores data in 
matrix so a solution can be computed. 

SMPT3 - This routine calls a routine to define 
additions to “B” matrix and “E” vector, plus 
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directing the called routine on where to store 
this information for the “B” matrix and ”E 
vector (for range and range difference data 


only). 


NETRAP - This routine uses a Newton-Raphson method to 
calculate time bias adjustment and then improve 
the receiver initial position. 


CROUT - This routine calls an IMSL routine that finds 
the determinate of a matrix. 


LINV3F - This is an IMSL routine that calculates the 
determinate of a matrix. 


POSAVG - This routine finds the average of all positions 
up to current time line. 


SMPT6 - This routine calls a routine to define additions 
to “B” matrix and “E” vector plus directing the 
called routine_on where to store this __ 
information (for range data only, satellite 
biases) 


LADRSC - This routine stores data in matrix so solution 
so solution can be computed after the range 
solution matrix and satellite clocks have béen 
updated. 


BEDI. - This routine eliminates certain parameters in 
the “B” matrix. 


BMSTR - This routine stores data in a matrix soa 
solution can be computed. 


SMIPT7 - This routine calls a routine to define 
additions to. the “B” matrix and “E” vector 
plus directing the called_routine on where to 
store this information (for range difference 
data only, satellite clock). 


EDRDSC - This routine stores data in a matrix so that 
a solution can be computed after updating the 
range solution matrix and satellite bias. 


SMPTS8 - This routine calls a routine to define additions 
to “B” matrix and “E” vector plus directing the 
called routine where to store this information 
(for range and range difference data only). 


SeIPplT9 - This routine calls a routine to define additions 
to the “B” matrix and “E” vector plus directing 
the called routine on where to store this 
information (for common time range data only). 


WDERSCE - This routine stores data in a matrix so that 
a solution can be computed after updating the 
Sonn range solution matrix and satellite 
clocks. 


Ey This routine finds data needed for elimination 
and then calls a routine to do that elimination. 


SMPT10 - This routine calls a routine to define 
aeaimonms to the 1 mutrix and vector plus 
directing the called routine on where to store 
this information (for Doppler data used as bias 
range data). 


a 


BSRNG - 


NEWITR - 


This routine defines additions to the “B” 
matrix ang) 4B “vecror 


This routine resets data so that a new iteration 
can be done. 


BLOCK DATA - This routine was added at NPS to make this 


program compatible with the FORTRAN 77 
conipilers on the IB M-3033 mainframe 
computer system. It contains common 
statements and data statements that 

NPS’s compilers would not permit where 
they were placed in the original code of 
program SOLTOM. 
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APPENDIX E 
INPUT GUIDE AND SOLUTION TYPES FOR PROGRAM SOLTOM 


( by S. L. Meyerhoff, NSWC) 


INPUT GUIDE 
Var Type Description 
CARD ONE Format (2F10.0,615) 


(This line can be read format free) 

col 1 to 10 

TSS real he GPS time tag fordatavaustant 
of solution (seconds 


COLNT tor20 
TES real The GPS time-tag of last data point 
to use in solution (seconds) 


col 21 to 25 (right justified) 
NPN integer The number of time lines of data 
needed before first solution 


26 to 30 (right justified) ie 
NEB integer The number of time lines of data 
between solutions 


col 31 to 35 (right justified) . 

IN ve integer The maximum number of iterations 
col 36 to 40 (right justified) . 

NSAT integer The maximum number of trackers to 


process at once (range | to 4) 


col 41 to 45 (right justified) 
IG] integer = | if data source is Standard 9-track 
format ; 
= 2 if data source GESAR test file 
format oe: ‘ 
= 3 if data Sounce 1suN7 a) eae 


col 46 to SO (right justified) _ 
Sie integer = 1 for satellite position from NAV 


ata 

= 2 for satellite position from 
Su Earth fixed trajectory 
iles 


col 51 to 55 (night justified) . 
ey Gre integer = | solve for clock aging, otherwise 
not solved for 


col 56 to 60 (right justified) 
IC EOGK integer = (0) solve clock corrections are from 
NAV data message 
= | clock corrections are post fit 
data from cards 


CARD TWO (receiver position) format (110,7G10.0 
(This card can be read format frec) 


col 10 
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fPeAG integer = ( for cartesian coordinatcs 
= | for geodetic coordinates 


Comhito 20 . 
LA | real latitude or X coordinate for position 
(degrees or kilometers) 
col 21 to 30 
ble) G real longitude or Y coordinate for position 
(degrees or kilometcrs) 
col 31 to 40 
i'l real height or Z coordinate for position 
(kilometers) 
col 41 to 50 
TB real time bias for receiver clock 
(kilomcters) 
col 51 to 60 
TD real trme drift for receiver clock 


(krlometers/second) 


col 61 to 70 
TO real epoch for clock drift (given_and 
solution) If negative, then first 
used timc line 1s the epoch for the 
solution clock drift 


col 71 to 80 
2 GIe) 


real aging for receiver clock 
(Kilometcrs/second*second) 
CARD THREE (Reference Geoid) Format {2G10.0) 


(This card can be read format frce) 
col 1 to 10 
AG 


real semi-major axis of the reference geoid 
(kilometers) 
Colsilito 20 -.. 
ROBL real. Obliquity of the reference geoid 
CARD FOUR (Sigmas for Position) Format (7G10.0) 
(This card can be read format frec) 
Comeau: 10 
SIGLT real sigma of X or latitude coordinate 
(Kilometers) 
Coll to 20 . 
SiG. G rcal sigma of Y or longitude coordinate 
(kilometers 
Coleleto 30 | 
SIGHT rcal sigma of Z or height coordinate 
(kilometers) 
col 31 to 40 
STG(1) real sigma of refraction scaling constant 
col 41 to 50 . . 
SIG) real sigma on recciver or satcllite | time 
bias (kilometers) 
col 51 to 60 
SIG(3) rcal sigma ol_recelver or satellite | 


clock drift (kilometers/second) 
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col 61 to 70 
SIG(4) real sigma for aging term or satellite 2 
bias (kilometers/second* second) 
CARD FIVE (Sigmas for Position) Format GY 
(This card can be read format free) 
colelatonl® 
S1G(5) real sigma on satellite 2 drift 
(kilometers/second) 
col 11 to 20 
S1G(6) real sigma for satellite 3 bias 
(kilometers) 
col 21 to 30 . 
SIG(7) real sigma for satellite 3 drift 
(kilometers/second) 
col 31 to 40 
SIG(8) real sigma for satellite 4 bias 
(kilometers) 
col 41 to 50 
SIG) real sigma for satellite 4 drift 
(kilometers) 


SOLUIION ] 1a 


CARD SIX (Types of Solutions to do) Format (1615) 
(This card cannot be read format free) 


Coles 
NSOL(1) integer = 1 for range solution receiver clock 
= no 
CaO 
NSOL(2) integer = 1 for ee difference solution 
receiver cloc 
Se) ive 
Collis 
NSOL(3) integer = | for range and range difference 
solution receiver cloc 
= Ono 
col 20 . 
NSOL(4) integer = | for common time range solution 
receiver clock 
= Ono 
cal 25 , | 
ICAL OW integer = i for Newton-Raphson range solution 
= 0 no 
col 30 
NSOL(6) integer = | for range solution satellite 
clocks 
= 0 no 
col 35 
NSOL(7) integer = Itfor range dillerence solution 
satellite clock 
= 0 6 
col 40 


NSOL(8) integer = | for range and range difference 


solution satellite clock 


== alte 
col 45 
NSOL(9) integer = | for common time range solution, 
Satellite clock 
= Qno 
col 50 
NSOL(10) integer = | for doppler data used as bias 


range data 
= ine 


CARD SEVEN (Variance Bias Factor) Format (2G10.0) 
(This card can be read format free) 


col 1 to 10 
K2 real bias factor for variance on range 
ata 
coll to 20 
K4 real bias factor for variance on doppler 
data 
CARD EIGHT (Print Options) acca) 
(This card can be read format free) 
cal 5 . - | 
IPRT(1) integer = | print position solution (X,Y,Z) 
= Q do not print 
col 10 
Pek i(2) integer = | print solution sigmas and clock 
terms 
= 0 do not print 
col 15 
IPRT(3) integer = | print data and time tags 
ean OMKGCOMpD Lol am oyatelt 
col 20 
IPRT(4) integer = 1 print residuals, range, trop, 
data sigma 
= 0 do not print 
col 25 . . im 
IPRT(S) integer = | print satellite position 
= 0 do not print 
Colo) ’ 
IPRT(6) integer = | print residuals and all data 
corrections 
= 0 do not print 
col 35 . | ee 
IPRT(7) integer = | print solution position in 
lat, long, height 
= 0 do not print 
col 40 
IPRT(8) integer = | print real time solutions from 
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9-track 
= 0 do not print 


CARD NINE (Weather Data, Elev. Cutoff) Format (4G 10.0) 
(This card can be read format free) 


col | to 10 7 

PB real pressure in millibars 
col 11 to 20 ' 

Te real temperature in degrees Centigrade 
col 21 to 30 tees: 

HD real relative humidity in percent 
col 31 to 40 

BE real elevation angle cutoff (degrees) 
CARD TEN (Station Number Format (15) 


(This card can be read format free) 


col 1 to 5 —— 
ISTA integer station identification 


CARD ELEVEN (Tracker 1 Clock Parameters) Format (3G 10.0) 
col 1 to 10 
AO 


real clock bias of satellite tracked on 
tracker 
col 11 to 20 . . 
Al real clock drift of satellite tracked on 
tracker 
col 21 to 30 . 
A2 real clock ane of satellite tracked on 
tracker 


CARD TWELVE (Tracker 2 Clock Parameters) Format (3G10.0) 
his card can be read format free) 


COlulpto7 LO 
AO real clock bias of satellite tracked on 
tracker 
col 11 to 20 
Al real clock drift of satellite tracked on 
tracker 
col 21 to 30 
A2 real clock ao of satellite tracked on 
tracker 


CARD THIRTEEN (Tracker 3 Clock Parameters) Format (3G10.0) 
us card can be read format free) 


col fto 10 
AQ 


real clock bias of satellite tracked on 
tracker 
col 11 to 20 
ral real clock drift of satellite tracked on 
{rackKer. 
col 21 to 30 
A2 real clock ee of satellite tracked on 
tracker 
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CARD FOURTEEN (Tracker 4 Clock AQeIEIES Format(3G10.0) 


eolpleto. 10 
AO 


real 
Cole to 20 
Al real 
col 21 to 30 
A2 real 


(This card can be read 


Hiatt ree) 


clock bias of satellite tracked on 
tracker 4 


clock drift of satellite tracked on 
tracker 


clock aging of satellite tracked on 
tracker 
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APPENDIA F 


KALMIN2 SUBROUTINE DESCRIPTIONS AND NESTING 


CARDIO - . 


DER. = 


SIGAOYZ = 


MATRIX - 


LINV2E - 


VMULFF - 


VMULFM - 


IND AT == 


TIMEUP - 


PROMPT - 
GEIDAT - 


NEWRAN - 


DEDPDT - 


DEFBR - 
CRDAT - 


BSADR - 


CRDAT? - 


ION” - 


IONRD - 
EPHDAT - 


CLOCK - 


This routine reads mnput for program. 

This routine converts latitude, longrtude and 
height to X, Y,.Z in earth-lixedeeocentric 
cartestan coordinates. 

This routine converts stgmas in the geodetic 
coordimate system into Siena in the earth- 
frxed geocentric coordinate system. 


This routine performs matrix operations by 
calling the proper routine mm the IMSL library. 


This routine is in the IMSL and performs a 
Ma UlIX Ally Crslonm: 


This routme ts m the IMSL and performs a 
matrtx multrplication. 


This routine performs a transpose of 
multiplrcatron of matrices. 


‘This routine grves the initial values to the 
B” matrix (U in Kalman filter) 


This routine uses a UDU’ factorization method 
to update state matrix. 


This routrne ards m the UDU’ factorization. 


This routine reads the next set of satellite 
tracking data from CON9TR output file. 


This routrne reduces norse and multipath im the 
range data. 


This routine converts Doppler phase mto range 
differenced data. 


Thrs routine defines bias range. 


This routine defines corrections for receiver 
1as. 


This routine will correct data for recerver - 
clock bras and drift and also correct GPS time- 
tag for the end of week cross over, 


This routine converts data mto range residuals 
and corrects data for error sources. 


This routine corrects for ronospheric 
refraction. 


This routine corrects range difference data 


This routine defines the values from the 
ephemeris data. 


This routine adds satellite and clock error and 
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Sa rrOS - 


EOLeOS - 


ZANG - 


hes 
SMOT - 


TROP - 


ROBE: 


ESiIG - 


GDOP - 


LINV3F - 


KALFIL - 


INKAL - 


UDFACT - 


RhAl  - 


UPDATE. - 


DEKAL-- 


CONSTR - 


AY Zee 


PRNSOL - 
COME L.- 


converts it to range measurements in kilometers 
ane also perform the same for range difference 
ata. 


This routine will find the position of the 
satellite from the navigation data message or 
from the Celestial Earth-Fixed Trajectory 
(Precise Ephemeris). 


This routine finds the position of the 
satellite from the Precise Ephemeris data. 


This routme finds the zenith angle by doting 
the X, Y, Z positions of the satellite receiver 
osition Vector to the vector from the receiver 
o the satellite. 


This routine finds the residual for ranges. 


This routine determines an earth rotation 
correction. 


This routine calculates a tropospheric 
refraction correction. 


This routine finds a relatrvistic correction 
to the data. 


This routine calculates the variance for raw 
scudorange and range differenced data for the 
satellite being tracked. 


This routine tf turned on, writes a GDOP file. 


This routine is in the IMSL and performs a 
matrix mnversion 


This routine uses corrected data to do the 
Kalman filter. 


This routine sets up mitral covariance matrix 
mew format. 


This routine puts Apriori covariance matrix 
mito UDU’ factorized matrix. 


This routine calculates partial derivatives of 
the data equation for pscudorange data. 


This routine will update the Kalman solution 
for one data point. 


This routine calculates partial derivatives of 
the data equation for delta (Doppler) range 
differences. 


This routine will constrain the solution to. 
the surface of a sphere (height) and constrain 
velocity. 


Convers X, Y, Zim amcarth-lixed ceocentric 
coordinate system to geodetic coordinate system. 


This routine prints solutions. 


This routine finds the covariance matrix from 
the UDU’ factorized matrix in U. 


BEOGK DATA Mine Toutine was added at INI’S to handle 
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data statements and common Statements. eel 
had to be moved in the code or written, 
be compatible with the [BM-3033 FOR’ rRAN T7 


compilers. 


ee PROGRAM KALMN2 


PbEN eZ 
SIGS 17 
MATRIX | 
PSV 2E 
Y NEC iar 
VMULFM 
INDAT 
TIMEUP 
PROMPT 
GETLDAL 
DEW 
a PBR 
CRDATI 
BSADR. 
C 2 
NEWRAN, 
IONRD 
ION 
IONRD 
ERGOT 
CLOCK2 
CLOCK 
SATO om 
> one O'S 
ZANG 
RES | 
SMOT 
ROP 
REL. 
FSIG 
GDOP ae 
LIN oie 
KAI 
INKAL 
J er 
LLIVa2 
KAI J 
UPDATE 
DPKAL 
UPDaAt i: 
CONSTR 
IPDATE 
or een 
PRN Ole 
COVEVL 


BLOCK DATA 
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APPENDIX G 


INPUT GUIDE FOR PROGRAM KALMN2 


Var Type 
CARD ONE 


( by S. L. Meyerhoff, NSWC) 
Description 
Format (2F10.0 


615) 


(This card can be read format free) 


COM to 10 
TSS 


real 
CGlmml to: 20 
TES real 
col 21 to 25 (right justified) 
NPN integer 
col 26 to 30 (right uses 
NPB Interer 
Com to 55 (right justied) 
ERY integer 
col 36 to 40 (right justified) 
NSAT integer 


col 41 to 45 (right justified) 
IGT integer 


col 46 to 50 (right justified) 
ISP integer 


col 51 to 55 (right justified) 
NEEOCK integer 


col 56 to 60 (right justified) 
IAGE integer 


col ue to 65 (right justified) 
Wa Ga) a integer 


CARD TWO (Receiver Position) 


The GPS time-tag for data at start 
of solution (seconds) 


The GPS time-tag of last data point 
to use in solution (seconds) 


The number of time lines between 
solutions 


The number of solutions between 
solutions being printed 


The number of solutions printed 
covariance matrix prints 


The number oi: ENG SOs uo i 
in the solution (range lI 


lec seurce 1s Standard 
ack format 
dais source 1s other 


vel 


= | for satellite Boe from 
NAVIGATION DATA MESSAGE 
= 72 for satellite. pesiien from 
Celestial Earth 

Trajectory files 


() clock corrections are from 
NAVIGA TION DATA M ESSAGE 
= | clock corrections are post 
fit data ipo cards 


= | make GDOP file 
Barre no file 1s made 


= | solve for clock aging 
otherwise no file 1s made 


Format (110,7G 10.0) 


(This card can be read format free) 


col 10 
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IFLAG integer = 0 for Cartesian coordinates 
= | for Geodetic coordinates 
col 11 to 20 
XLAT real latitude or X coordinates for 
position (degrees or kilometers) 
col Jie 50 
VY EOG real longitude or Y coordinate for 
position (degrees or kilometers) 
col 31 to 40 
Zann real height or Z coordinate for 


position (kilometers) 
col 41 to 50 
TB 


real time bias for receiver clock 
(kilometers 
col SI to 60 . 
iD real time drift for receiver clock 


(kilometers/second 


col 61 to 70 
TO real epoch for clock drift (given. and 
solution) if negative, then first 
used time line 1s the epoch for 
the solution clock drift 


col 71 to 80 
AGE aging for receiver clock 
(Kilometers/second* second) 


CARD THREE (Reference Geoid) Format (2G 10.0) 
(This card can be read format free) 


real 


col I to 10 aa 
AG real semi-major axis of the reference 
geoid (kilometers) 


COME ton) ae 
ROBL real obliquity of the reference geoid 


CARD FOUR (Sigmas for Position) Format (7G10.0) 
(This card can be read format free) 


col 1 to IO ; | 
SIGES real sigma of X or latitude coordinate 
(kilometer) 
Colmisiance 
SIGLG real sigma of Y or longitude coordinate 
(kilometer) 
col 21 to 30 
SIGE real sigma of Z or height coordinate 
(kilometers) 
col 31 to 40 
STG(1) real sigma on refraction scaling constant 
col 41 to 50 
STG(2) real sigma on receiver time bias 
(kilometer) 
col 5I to 60 
SIG(3) real sigma on receiver clock drift 


(kilometer/secon 
col 61 to 70 


70 


S1G(4) 


CARD FIVE (Sigmas for Velocity) ) 
(This card can be read format free) 


col 41 to 50 
SIG) 


CArse SIX (Types af solute 


real 


real 


real 


real 


real 


real 


integer 


integer 


integer 


integer 


integer 


integer 


integer 


integer 


sigma of velocity in X direction 
(Kilometers/second) 


format (7G10.0 


sigma of velocity in Y direction 
(kilomcters/second) 


sigma of velocity in Z direction 
(kilometers/second) 


sigma On receiver time bias 
(kilometers) 


sigma on receiver clock drift 
(kilomcters/second) 


sigma not used 


ns todo) Format (1615 
his card cannot be read format free) 


= | for range position solution 
= 0 no 


= | Height constraint on position 
= () constraint only if less than 4 
satellites 


] write position solution on file 
0 no solution save 


= | to solve for refraction 
correction in position 
solution 

= Qno 


= | to solve for refraction 
correction in velocity 


solution 
= Qno 


= | to solve for clock drift 
from position solution 
OQ no 


= | doa velocity solution 
this has not been used 
efore) 

= Ono 


= | doa time update of 
covariance matrix 
me ne 
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col 45 
NSOL(9) integer 


col 50 
NSOL(19) integer = | for a constraint on 
velocity solution 


= 0 no 


CARD SEVEN (Variance Bias Factor) Format (2G10.0) 
(This card can be read format free) 


col 1 to 10 
K2 real bias factor for variance on range 
data 
col dl to.z0 
K4 real bias factor for variance on 


Doppler data 


CARD EIGHT (Print Options) Format (515) 
(This card cannot be read format free) 


col 5 
IPRT(I integer = I print position solution 
(1) g xv Pp | 
= 0 do not print 
col 10 | . —: 
eZ) integer = | print solution sigmas and 
clock terms 
= Q do not print 
Colmes 
PPT (3) integer = | print data and time tags 
= 0 do not prm 
col 20 
IPRT(4) integer = | residuals, range, trop, 
data sigma ss 
= 0 do not print 
col 25 , ms 
IPRT(S) integer = | print satellite position 
= 0 do not print 
col 30 
IPRT(6) integer = | residuals and all data 
corrections 
= Q do not print 
col 35 | ee 
Pe Rae) integer = | print position in lat, long, 
and height 
col 40 | | | | 
IPRT(8) integer = | pa real time solutions from 
-track 


CARD NINE (weather Data, Elev. Cutoff) Format (4G10.0) 
(This card can be read format free) 


col | to 10 7 
PB real pressure in millibars 


col 11 to 20 
PC real temperature in degrees Centigrade 


eZ 


Gori to 30 
HID real relative humidity in percent 


col 31 to 40 
El real elevation angle cutoff (degrees) 


CARD TEN (Station Number) Format (Free) 
(This card can be read format free) 


col I to 10 
QI) 


real time pro aloe noise for X 
velocity bk ometers) 
col I1 to 20 
Q(2) real time propagation noise for Y 
velocity (kilometers) 
Cole to 30 
(3) real time propagation noise for Z 
velocity (kilometers) 
col 31 to 40 ; ; ; 
Q(4) real time propagation noise for time 
drift 
col 41 to 50 _ 
SIGPS real sigma On position constraint 
col 51 to 60 . 
SIGYL real sigma on Velocity constraint 


CARD ELEVEN (Diagonal elements of state transition matrix) 
1s card can be read format free) 


Coimaleto 10 7 
PHI(1I,1) real state transition matrix for 
parameter 
col 1] to 20 _ 
Pivl(2,2) real State transition matrix for 
parameter 2 
col 2] to 30 . 
Feel 3,3) real state transition matrix for 
parameter 
col 31 to 40 _ 
PH1(4,4) real state transition matrix for 
parameter 4 
col 41 to 50 - . 
Peis ;)) real State transition matrix for 
parameter 5 
col 51 to 60 a2 
PH1(6,6) real State transition matrix for 
parameter 
coleo! to. 70 a 
Fit 7,7 real state transition matrix for 


parameter 7 
CARD TWELVE (Diagonal elements of state transition matrix) 
us card can be read format {ree) 
ol s . 
PIII(8,8) real state transition matrix for 
parameter 8 


Col ieton 2) 


13 


Feri 9,9) real state transition matrix for 
parameter 9 


CARD THIRTEEN (lleight HOF position constraint 
(This card can be read format free) 


col | to 10 
rie real Height above reference 
ellipsoid to constrain 
osition solution to 
kilometer) 
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APPENDIX H 
ABBREVIATIONS AND ACRONYMS 


Giles = Global! Positioning System 

NPS = Naval Postgraduate School Monterey CA 

DMA = Defense Mapping Agency 

NOAA = National Oceanic and Atmospheric Administration 

NSWC = Naval Surface Weapons Center 

SV = Space Vehicle 

wo = Spherical Error Probability 

P-code = Precision code 

C/A-code = Coarse Acquisition code 

WGS = World Geodetic System 

PDOP = Position Dilution of Precision 

HDOP = Horizontal Dilution of Precision 

VDOP = Vertical Dilution of Precision 

TDOP = Time Dilution of Precision 

Sie = Standardized Exchange Format 

CIO = Conventional International Origin 

ARL = Applied Research Laboratories at the University of Texas 
at Austin 

IBM = International Business Machines 

CDEC = Computer Data Corporation 

m = meters 

kin = kilometers 

S = second 

EZ = Hertz 

Miz = Megahertz 
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